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PREFACE 


A summary of equations often used in free -flight and wind-tunnel data reduction and 
analysis is presented. Included are transfer equations for accelerometer, rate -gyro, and 
angle -of -attack instrumentation; axes -system transfers of aerodynamic derivatives; and 
methods for measuring moments of inertia. In general, the equations are in a complete 
form; for example, those terms are retained that are missing when planar symmetry is 
assumed for airplanes. 
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INTRODUCTION 


The equations in this report are the coordinate transformation and motion equations 
used in the various tasks associated with free -flight and wind-tunnel data reduction and 
analysis. These tasks range from reducing flight data to calculating the motions on a * 
digital or analog computer and to applying various techniques for analyzing the data, such 
as in references 1 and 2. 

While many publications contain a number of these equations, no one contains all 
that are usually needed in a complicated aerodynamic analysis; even the more nearly com- 
plete reports (refs. 3 and 4, for example) omit the equations for transferring aerodynamic 
stability derivatives from one moment reference to another. Moreover, in most cases 
the equations are simplified, when they are presented, by assumptions such as small 
angles of attack, zero sideslip, and small perturbation motions. Expanded forms of many 
of the equations, on the other hand, are needed in special problems that may arise. For 
example, parawing vehicles, which have their center of gravity located well below the wing 
surface, require the expanded forms of the axes transformations when data measured about 
a point on the wing are to be transferred to the center of gravity; reentry motion studies 
sometimes involve large -amplitude motions so the complete forms of the transformations, 
without the assumptions of small angles of attack or sideslip, are needed. The engineer 
working on any of these special problems usually has to derive these equations himself, 
and this can be time consuming. 

The purpose of this report is to provide the basic equations from which many of the 
equations needed in a particular analysis can be generated. A comprehensive summary 
of the basic axes transformation and motion equations is included, with most of these 
given in their expanded, most general forms. Once these expanded forms are available, 
the simpler forms can be written out fairly easily, and yet the general forms are here 
when needed for special cases. 

The general forms presented include axes transformations that enable transfer back 
and forth between any of the five axes systems that are encountered in aerodynamic analy- 
sis. Equations of motion are presented that enable calculation of motions anywhere in the 
vicinity of the earth. Special problems are also considered; since flight instruments, 
such as accelerometers or rate gyros, are not always alined along mutually perpendicular 
axes, the procedure for correcting instrument readings for nonorthogonal alinements is 
outlined. 

In addition to these general forms, many of the simplified forms used frequently in 
practical applications are summarized in appendix A. 

Other relationships are presented in appendixes B to F. For example, appendix C 
summarizes methods used to measure moments of inertia of models and full-scale 



vehicles; appendix E discusses the use of the direction-cosine and the quaternion methods, 
often used in place of Euler angles in specifying vehicle alinement; appendix F discusses 
the scaling parameters used in model testing. However, throughout this paper, the empha- 
sis is on providing the basic equations. For discussions of their development and of the 
procedures used in their application, the user should turn to general published works on 
flight-motion analysis. A comprehensive bibliography of these works is provided and 
includes textbooks and reports dealing with stability, control, and performance as well 
as reports discussing various techniques for extracting stability derivatives from flight 
data. 

SYMBOLS 


Throughout this paper, symbols are defined in terms of SI Units with equivalent U.S. 
Customary Units given parenthetically. Factors for converting from U.S. Customary to 
SI Units are given in table I. 

{ generalized angle of attack, defined for various axes systems in tables III 
and IV, rad (deg) 

cross-sectional area in eq. (F-ll) of appendix F, m^ (ft^) 

A f angle between surface wind vector and plane of local horizontal, measured 

perpendicular to plane of local horizontal (fig. 8), rad (deg) 

Aj) acceleration along flight path, g units (lg = 9.807 m/sec^) 

A^ acceleration in lift direction, g units (lg = 9.807 m/sec2) 

A v ,A V o^^Afv components of acceleration along X,Y,Z vehicle reference axes 

A, Cg’ 1 ,Cg’ Zj,Cg 

at c.g., respectively, g units (lg = 9.807 m/see*) 

A x i? A Y i,A z i components of acceleration indicated by accelerometers along 

X^,Y^,Z^ instrument axes, respectively, g units 
(lg - 9.807 m/sec^) 

a speed of sound, m/sec (ft/sec) 

a r damping constant defined by eq. (D-4) of appendix D, dimensionless 

a diameter of circle around which wires or rods are attached in bifilar or 

trifilar methods of measuring moments of inertia in eq. (C-9) of 
appendix C (see, also, fig. 14(b)), m (in.) 
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a e 

a s 

B 

B* 

b 

C 


equatorial radius of earth’s reference ellipsoid, m (ft) 

distance from knife edge to spring (fig. 13), m (ft) 

generalized angle of sideslip, defined for various axes systems in tables IIP 
and IV, rad (deg) 

angle between projection of wind vector onto plane of local horizontal and 
component of wind velocity tangent to circle of local horizontal (fig. 8), 
rad (deg) 

span, m (ft) 

constant in eq. (D-3) of appendix D, rad (deg) 


C A 

C C 

C D 

C T 

Cl 


axial -force coefficient 
crosswind coefficient 
drag coefficient 

drag coefficient for wind-tunnel stability axes 
lift coefficient 


lift coefficient in steady -state level flight 


c*,c m ,c n 


generalized aerodynamic moment coefficients about X,Y,Z vehicle 

„ . Moment 

reference axes, respectively, — - — — — 

body -axes moment coefficients 


Cm,o 

C N 


pitching -moment coefficient in steady -state level flight 
normal -force coefficient 


C X> C Y’ C Z 


generalized aerodynamic force coefficients about X,Y,Z vehicle refer 

Force 

ence axes, respectively, — - 

q S 

A oO 

body -axes force coefficients 
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c reference chord, m (ft) 

Cp specific heat at constant pressure, j/kg-°C (Btu(thermochemical)/lbm-°F) 

c v specific heat at constant volume, j/kg-°C (Btu(thermochemical)/lbm-°F) 

/aerodynamic drag in eqs. (1-25) of Sec. I, N (lbf) 

(direction -cosine matrix defined in eq. (E-2) of appendix E 

d X> d Y’ d Z displacements of centroid of jet interface, measured with respect to X,Y,Z 
vehicle reference axes, respectively, m (ft) 

E modulus of elasticity, N/m 2 (lbf/ft 2 ) 


eQ,e-|.,e2,eg Euler parameters in eq. (E-9) of appendix E 

F r,X’ F r,Y’ F r,Z forces due to jet control (not main rockets) along X,Y,Z vehicle 

reference axes, respectively, N (lbf) 

F x ,F y > F Z forces along X,Y,Z vehicle reference axes, respectively, N (lbf) 

G transformation matrix composed of Euler parameters 

GE geocentric gravitational constant, m 3 /sec 2 (ft 3 /sec 2 ) 

g acceleration due to gravity, m/sec 2 (ft/sec 2 ) 


^X’^Y’^Z components of gravitational acceleration along Xg,Yg,Zg gravity axes, 
respectively, for an oblate earth, m/sec 2 (ft/sec2) 


H angular direction vehicle is traveling (fig. 6); angle measured positive 

clockwise from north, i.e., for H = 90°, vehicle is traveling toward east, 
rad (deg) 

h altitude, m (ft) 

I moment of inertia, kg-m 2 (slug-ft 2 ) 
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o 

I mounting -cradle moment of inertia about vehicle roll axis, kg-m 

(slug-ft2) 

I rm moment of inertia of engine rotating mass, kg-m 2 (slug -ft 2 ) 

I X ,I Y ,I Z moments of inertia about X,Y,Z vehicle reference axes, respectively, 
kg-m 2 (slug-ft2) 

I XY’ I XZ’ I YZ products of inertia, 1^ = J xy dm, I xz = J xz dm, 

*YZ = yz dm, kg-m 2 (slug -ft 2) 

i,j,k unit vectors having properties given in eqs. (E-10) of appendix E 

J 2 ,J 3 ,J 4 second, third, fourth zonal harmonics, respectively, of earth’s reference 
ellipsoid (values given following eq. (V-45) of Sec. V) 

j unit vector (see i,j,k) 

K spring constant, N/m (lbf/in.) 

Ki,K 2 ,...,Kg constants defined in eqs. (1-13) of Sec. I 

K 0 boom misalinement angle at zero g, measured with respect to X vehicle 

reference axis; positive with boom inclined above X-axis, rad (deg) 

K s constant for springs shown in fig. 13, N/m (lbf/ft) 

K g i,K g 2 values of K s for springs of different strengths, N/m (lbf/ft) 

Kj. 1 ,K t 2 constants, determined from K s i and K g 2 , used in measuring moments of 
inertia in eqs. (C-6) and (C-7) of appendix C, N/m (lbf/ft) 

unit vector (see i,j,k) 

torsion-spring constant in eq. (C-8) of appendix C, N -m/rad (ft-lbf/rad) 

k ( 

coefficient of heat conduction in eq. (F-9) of appendix F, J/m-sec-K 
(Btu(thermochemical)-in./ft2-sec-°F) 
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L,M,N 


(aerodynamic lift in eqs. (1-25) of Sec. I, N (lbf) 

\ geocentric latitude of vehicle, positive in northern hemisphere, in Sec. V 
and appendix E, rad (deg) 

aerodynamic moments about X,Y,Z vehicle reference axes, respectively, 
N-m (ft -lbf) 


L M N 
^rm’Tm^rm 


moments due to gyroscopic action of engine rotating mass about 
X,Y,Z vehicle reference axes, respectively, N-m (ft -lbf) 


ft' 


/length used in nondimensionalizing moments, m (ft) 


/ 


[characteristic length, m (ft) 

length of wire or rod on which model is suspended in measuring moment of 
inertia by multifilar -pendulum method in appendix C, m (ft) 




'"distance between fore and aft weighing scales in eq. (C-l) of appendix C, 
m (ft) 


distance from knife edge to supporting spring in eq. (C-5) of appendix C, 
_ m (ft) 


^XX’ £ XY’-*-’ £ ZZ 


direction cosines defining orientation of rotated (primed) axes 
system with respect to initial (unprimed) axes system; &xx * s 
cosine of angle between X-axes of two systems; #XY> cos i n e of 
angle between X-axis of rotated system and Y-axis of initial sys- 
tem; . . .; and J^zZ’ cosine of angle between Z-axes of two sys- 
tems (eqs. (E-3) of appendix E) 


M Mach number 

M aerodynamic moment (see L,M,n) 

M rm moment due to engine rotating mass (see L r m> M rm?Nrm) 

M r X’^r Y’-^r Z moments due to jet controls (not main rockets) about 

X,Y,Z vehicle reference axes, respectively, N-m (ft-lbf) 


M^,My,M 2 moments about X,Y,Z vehicle reference axes, respectively, 

N-m (ft-lbf) 


6 



m 


instantaneous mass of vehicle, kg (slugs) 


m’ = 


mVoo 


m 


m c 

N 

Nrm 

P 


Pi 

p 2 

P M 

P M+S 

P S 


Pq^o^o 


P 


P,q,r 


Pi J< ^i’ r i 


Qo 


time rate of change of vehicle mass (such as due to fuel consumption); 
negative when vehicle losing mass, kg/sec (slugs/sec) 

mass of weighing cradle used in measuring moments of inertia, kg (slugs) 

aerodynamic moment (see L,M,n) 

moment due to engine rotating mass (see Lrm>M rm ,N rm ) 
period of oscillation, sec 

period of oscillation with spring constant Kj. sec 
period of oscillation with spring constant Kj- 2> sec 
period of oscillation for model alone, sec 

period of oscillation for model plus supporting hardware, sec 

period of oscillation for supporting hardware, sec 

initial, steady -state values of angular -velocity components along 
X,Y,Z vehicle reference axes, respectively, rad/sec 

static pressure, N/m^ (lbf/ft2) 

components of angular velocity about X,Y,Z vehicle reference axes, 
respectively, rad/sec (deg/sec) 

components of angular velocity indicated by rate gyros about 
Xj,Y|,Zj_ instrument axes, respectively, rad/sec 

angular -velocity component (see Po?Qo> p o) 



angular -velocity component (see p,q,r) 
quaternion in eq. (E-9) of appendix E 


impact pressure, N/m 2 (lbf/ft 2 ) 
angular -velocity component (see Pi,qi,r^ 
free -stream dynamic pressure, N/m^ (lbf/ft 2 ) 
conjugate of quaternion in eq. (E-12) of appendix E 

[distance of c.g. of vehicle from center of earth (fig. 4), m (ft) 

(scale factors, with subscripts, as defined in table VI and discussed 
in appendix F 

radius of assumed spherical earth, 6378.123 km (20 925 631 ft) 
angular -velocity component (see P 0 ,Q 0 ,R 0 ) 

vertical reaction force (fig. 10) in eq. (C-2) of appendix C, N (lbf) 

forces at supports in setups for determining c.g. location in eq. (C-l) of 
appendix C, N (lbf) 

angular -velocity component (see p,q,r) 

vector expressing distance from c.g. to origin of instrument -axes system, 
m (ft) 

distance from origin of gravity -axes system to vehicle reference center, 
m (ft) 

angular -velocity component (see p^,q^,r^ 
vehicle reference area, m 2 (ft 2 ) 


{ total vehicle thrust in eq. (V-36) of Sec. V and in eq. (A-134) of appendix A, 
N (lbf) 

temperature in eq. (F-9) of appendix F, K (°F) 



component of thrust along flight path, N (lbf) 


/Initial thrust in eq. (V-36) of Sec. V 

° ^representative temperature of flow in eq. (F -10) of appendix F, K (°F) 

T T T Z components of thrust along X,Y,Z vehicle reference axes, respectively, 
N (lbf) 


t time, sec 

t^^ time to damp to one -half amplitude, sec 

t-^ arbitrary time defined in fig. 15, sec 

t2 arbitrary time defined in fig. 15, sec 

U 0 ,V 0 ,Wo initial, steady -state values of linear -velocity components along 

X,Y,Z vehicle reference axes, respectively, m/sec (ft/sec) 

u,v,w components of vehicle absolute (inertial) velocity along X,Y,Z vehicle 

reference axes, respectively, m/sec (ft/sec) 

u,v,w components of relative velocity (velocity of vehicle with respect to air), 

along X,Y,Z vehicle reference axes, respectively, m/sec (ft/sec) 

u^v^w^ indicated components of relative velocity along X^,Y^,Z^ instrument axes, 
respectively, m/sec (ft/sec) 

V resultant velocity, m/sec (ft/sec) 

AV increment in resultant free -stream velocity, m/sec (ft/sec) 

V 0 velocity component (see Uo?V 0 ,W Q ) 

free -stream reference velocity, m/sec (ft/sec) 

v x >v Y ,Vr ^ components of geostrophic-wind velocity, due to earth rotation 

(atmosphere rotates with earth), relative to X,Y,Z vehicle reference 
axes, respectively, m/sec (ft /sec) 

v inertial -velocity component (see u,v,w) 

v relative -velocity component (see u,v,w) 
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Vi 

W 

W 

W c 


W M 


W M+S 

W S 

W 0 

w x ,w Y ,w 


w 

w 


w i 

X,Y,Z 

X,Y,Z 


x e 

x g , Y g, Z g 


x,y,z 


relative -velocity component (see Ui,Vi,WjJ 
vehicle weight, N (lbf) 

resultant velocity of earth surface wind, m/sec (ft/sec) 
weight of weighing cradle, N (lbf) 
weight of model alone, N (lbf) 

weight of model plus supporting hardware, N (lbf) 
weight of model supporting hardware, N (lbf) 
velocity component (see U 0 ,V o ,W 0 ) 

, components of wind velocity, due to surface winds, relative to 
X,Y,Z vehicle reference axes, respectively, m/sec (ft/sec) 

inertial -velocity component (see u,v,w) 

relative -velocity component (see u,v,w) 

relative -velocity component (see Ui,Vi,Wj() 

vehicle reference axes 

aerodynamic forces along X,Y,Z vehicle reference axes, respectively, 

N (lbf) 

right-handed inertial axes with origin at center of earth (fig. 4) 
gravity axes with origin at surface of earth (fig. 4) 

axes in orthogonal coordinate system with origin at an instrument at point 
x,y,z and with system alined at angles ^,6,(p with X,Y,Z vehicle 
reference axes at c.g., respectively, as shown in fig. 1 

distances measured from X,Y,Z vehicle reference axes, respectively, m 


t 

\ 


(ft) 
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x,y,z 

x g ,y g ,z g 

x X’^X’ z X 

x Y? y Y » z Y 

X Z’^Z’ Z Z 




J M+S 


distances from c.g. to another point on vehicle, m (ft) 

displacements of vehicle with respect to X g ,Y g ,Z g gravity axes, respec- 
tively (figs. 4 and 5), m (ft) 

* 

coordinates of X-axis accelerometer or rate gyro, measured from 

X,Y,Z vehicle reference axes at c.g., respectively , m (ft) 

coordinates of Y-axis accelerometer or rate gyro, measured from 
X Y Z vehicle reference axes at c.g., respectively, m (ft) 

coordinates of Z-axis accelerometer or rate gyro, measured from 
X,Y,Z vehicle reference axes at c.g., respectively, m (ft) 

distance from knife edge to vertical reaction force R x (fig. 10), m (ft) 

distance from knife edge to vehicle c.g. (fig. 13), m (ft) 

vertical distance to c.g. of weighing cradle (figs. 10 and 13), m (ft) 

distance from knife edge to c.g. of model, m (ft) 

distance from knife edge to c.g. of model plus supporting hardware, m (ft) 
distance from knife edge to c.g. of supporting hardware, m (ft) 


a 


AO! 


Aa. 


Aot< 


angle of attack, rad (deg) 

i w b 

tan" 1 — = tan 


w b W b ~ V Z,b ~ W Z,b 


u b 


u b " V X,b " W X,b 


position error, error due to location of angle -of -attack sensor in flow of 
vehicle, equal to free -stream angle of attack minus true or corrected local 
angle of attack; positive for a free -stream angle of attack greater than 
local angle of attack, rad 

difference between angle of attack and trim angle of attack at time t v rad 
difference between angle of attack and trim angle of attack at time t2 , rad 
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increment in angle of attack caused by boom bending under static air load, 
measured with respect to boom center line; positive for upward deflection 
of boom, rad 


change in angle of attack caused by bending of mounting boom under inertial 
load, measured with respect to boom center line; positive for upward 
deflection of boom, rad 

change in angle of attack caused by vehicle rotation, measured with respect 
to vehicle reference axis, rad 


angle of attack indicated by vane or other sensor, measured with respect to 
center line of mounting boom; positive with nose of vane pointed down, rad 


angle between X body axis and X flight stability axis (fig. 3), rad (deg) 
trim angle of attack, rad (deg) 

change in angle of attack caused by upwash from mounting boom* measured 
with respect to free -stream velocity; positive for upwash (upward flow 
normal to free stream), rad 

vane floating angle, angle caused by slight mass unbalance or by warp in 
vane surface due to imperfections in manufacture, measured with respect 
to local velocity vector; positive for nose -up deflections of vane, rad 


angle of sideslip, rad (deg) 

!] 

V 


sin-1 lb = gln -l V b - V Y,b - W Y,b 




sideslip position error, rad 


increment in sideslip angle caused by boom bending under static air load, 

measured with respect to boom center line; positive for boom deflected 
to right, rad 

change in sideslip angle caused by bending of mounting boom under inertial 
loads, measured with respect to boom center line; positive for boom 
deflections to right, rad 



3 d change in sideslip angle caused by vehicle rotation, measured with respect 

to vehicle reference axis, rad 

s . sideslip angle indicated by vane or other sensor, measured with respect to 

center line of mounting boom; positive with nose of vane pointed to right ^ 

of boom, rad 

change in sideslip angle caused by sidewash from mounting boom, measured 
with respect to free -stream velocity; positive for positive sidewash, rad 

0 sideslip -vane floating angle, angle caused by warp in sideslip vane due to 

imperfections in manufacture, measured with respect to local velocity 
vector; positive for vane deflected to right, rad 

T transformation matrix for orthogonal axes system, defined in eq. (A-2) of 

appendix A 

■^nonorthogonal transformation matrix for nonorthogonal axes system, defined in 

eq. (1-14) of Sec. I 

(ratio of specific heats in eqs. (1-19) and (1-20) of Sec. I, c p /c v 

^ ^flight -path angle in eq. (V-17) of Sec. V (fig. 6), rad (deg) 

6 control deflection, rad (deg) 

5 5 5 control deflections (aileron, elevator, rudder, respectively) rad (deg) 

a> e> r 


5rPM change in engine rpm 

e angle between Xp principal and X body axes (fig. 3), rad (deg) 

strain (elongation per unit length) measured on model, cm/cm (in. /in.) 

6 strain (elongation per unit length) measured on prototype, cm/cm (in. /in.) 

P 

£ angle between X s flight stability and Xwt wind-tunnel stability axes 

(fig. 3), rad (deg) 

7] angle between Xp principal and Xg flight stability axes (fig. 3), rad (deg) 
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geocentric longitude of vehicle, positive counterclockwise looking in direction 
of positive Z e inertial axis (fig. 4) , rad 


li coefficient of viscosity, N-sec/m 2 (slugs/ft-sec) 

v kinematic viscosity in eq. (F-9) of appendix F, p/p, m 2 /sec (ft 2 /sec) 

£ phase angle, rad (deg) 

P atmospheric density, kg/m 3 (slugs/ft 3 ) 


f range angle in eq. (V-21) of Sec. V, rad (deg) 

(surface tension in eq. (F-7) of appendix F, N/m (lbf/ft) 


0 roll angle defined in fig. 10, rad (deg) 

0,0,0 Euler angles defining angular alinement of one axes system with respect 
to another axes system, rad (deg) 

^g’^g’^g Euler angles defining alinement of X,Y,Z vehicle reference axes, respec- 
tively, with respect to gravity-axes system (fig. 5), rad (deg) 


^rm> 0 rm Euler angles describing alinement of engine -thrust axes (fig. 9), rad (deg) 

^x^x ,0 X Euler angles of X-axis accelerometer or rate gyro, measured with respect 
to X,Y,Z vehicle reference axes at c.g., respectively, rad (deg) 


0Y?0Y>0y Euler angles of Y-axis accelerometer or rate gyro, measured with respect 
to X,Y,Z vehicle reference axes at c.g., respectively, rad (deg) 

^Z’^Z’^Z angles of Z-axis accelerometer or rate gyro, measured with respect 

to X,Y,Z vehicle reference axes at c.g., respectively, rad (deg) 


rate of notation of earth, rad/sec 


fi rm 


angular velocity of engine rotating mass, rad/sec 


co 
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resultant angular velocity in Sec. I, rad/sec 
frequency of oscillation in appendix D, rad/sec 



Ci) X t >a) Y r >U) Z T angular -velocity components along X T ,Y T ,Z f vehicle reference axes, 
respectively, usually measured with respect to gravity -axes or 
inertial-axes system, rad/sec 

Subscripts: • 

b body axes 

o initial conditions 

P principal axes 

s flight stability axes 

w wind axes 

wt wind-tunnel stability axes 

Notation: 


. (Dot) first derivative with respect to time 

1 (Prime) unless otherwise specified, a primed quantity is one referred to axes sys- 
tem located at point x,y,z and/or alined at angles with respect 

to initial reference axes system 

A perturbation quantity (unless specified otherwise) 

Subscripts used with coefficient symbols denote derivatives as follows: 


oi 

a 

0 

4 

u 


with respect to a 
with respect to a H/2Voo 
with respect to fi 
with respect to 
with respect to u/Voo 
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with respect to v/Voo 


with respect to 
with respect to 
with respect to 
with respect to 


w/Vso 
p l/2V x 
q^/av*, 

r£/2V«, 


with respect to control deflection 

with respect to V/Voo. This derivative, for a force such as the X-force, is 

9C 

equal to 2Cv + - , where the term 2C X represents the change in 

A v/Voo) ac x 

X-force due to changes in free -stream velocity and the term 

represents the change in C X due to effects such as Mach number or 
aeroelastic effects. If these effects are negligible, then C X y = 2C X , 
Cmy = ^ Cj-q , etc. 



SECTION I 


EQUATIONS INVOLVING BASIC FLIGHT MEASUREMENTS 


This section summarizes some of the relationships used in working with basic flight 
measurements. Included are general equations for transferring accelerations, linear 
velocities, and angular velocities between any two axes systems on a flight vehicle that 
are related by a 10,0,0 Euler rotation sequence; the simplified forms that are frequently 
used to correct accelerometer and rate -gyro readings for instrument displacement and 
misalinement are given as equations (A-5) to (A-16) in appendix A. Also included are 
the correction equations for angles of attack and sideslip, equations for determining flight 
Mach number from measurements of static and impact pressures, and equations for deter- 
mining vehicle forces and moments from accelerometer and rate -gyro readings. 

GENERAL AXES TRANSFORMATIONS FOR COMPONENTS 
OF ACCELERATION, LINEAR VELOCITY, AND 
ANGULAR VELOCITY 

The general equations given here transfer between the two axes systems shown in 
figure 1, where the "instrument -axes" system represents any system displaced and/or 
misalined with respect to the vehicle reference axes at the c.g. It should be noted that the 
^i>^i?^i components of the instrument -axes system need not be mutually perpendicular 
(orthogonal) nor referred to the same origin. This would be the case, for example, if 
accelerometers were used to measure the separate components of acceleration but the 
accelerometers were not orthogonally alined and were located at different points on the 
vehicle. Hence, three separate instrument axes have to be considered. For the equa- 
tions developed here, one instrument is assumed to be alined along the X-axis of a sys- 
tem located at point ^X’^X’^X an< ^ alined at Euler angles ^X’^X’^X w h:h the vehicle 
reference axes; a second instrument is assumed to be alined along the Y-axis of a sys- 
tem located at point anc * a * in ed at Euler angles tyy,9y,<fiy w ^h the refer- 

ence axes, and so on. Nine displacement coordinates and nine Euler angles in all are 
needed to define the locations and alinements of the three instruments. 

Both axes rotation and axes translation can be performed by using these 

equations ^in the acceleration equations, axes -translation terms are of the form 
— dV — 

a inertial “ ~[j“ + w x V; in the velocity equations, they are of the form 

^inertial ~ + w x r j. The transformations from the vehicle reference axes to the 

instrument axes are given in equation form; however, those from the instrument axes 
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SECTION I - Continued 


to the vehicle reference axes are given only in matrix form , since they involve a matrix 
inversion which becomes long and complicated unless certain simplifying assumptions 
are made. Two simplified cases, in which (1) instrument axes are orthogonal and (2) 
misalinement angles are small, are given as instrument -correction equations in equa- 
tions (A -5) to (A-16) of appendix A. 


Transfer From Vehicle Reference Axes to Instrument Axes 
Acceleration. - 


gA X,i 


gA x cg - (r 2 + q 2 )x x + (pq - r)y x + (rp + q)zJcos 0 X cos ^ x 
+ jgA y cg + (pq + r)xy - (p 2 + r 2 )y y + (qr - p)zyjcos 6 X sm ^ x 
- [ gA Z,cg + ( pr " q )*Z + ( qr + p )^Z “ ( q2 + p2 )*z] sin 6 X 


(I-D 


gA Y) i = [gA x>cg - (r 2 + q 2 )x x + (pq - r)y x + (rp + q)z x ](cos sin 0 y sin 0 y - sin * y cos 4> y ) 

+ jgA y cg + (pq + r)x y - (p 2 + r 2 )y y + (qr - p)z^J(sta ^ y sin 0 y sin 0 y + cos i// y cos 4> y ) 
+ [gA z cg + (pr - q)x z + (qr + p)y z - (q 2 + P 2 )z z ]cos 6 y sin 4> y 


(1-2) 


S A Z i = [g A x cg " ( r2 + q2 )^X + (pq " ^X + ^ rp + q ^X_ ( COS sin 9 Z cos <t>z + Sin sin .^z) 

+ jgA y cg + (pq + r)x y - (p 2 + r 2 )y y + (qr - p)z^J(sin sin 6 Z cos 0 Z - cos ^ z sin ^z) 

+ [g Az cg + (pr - q)x z + (qr + p)y z - (q 2 + p^z^cos 6% cos <p z (I" 3 ) 

Linear velocity . - 

u i = (u + qz x - ry x ^cos 6 X cos \p x + (v + rx y - pz y ^cos $ x sin \j/ x 
- (w + py z - qx^sin 6 X 
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+ qz x - ry x ^sin 0y sin C0S ^Y “ sin ^Y C0S ^y) 

+ + rXy " pZy^sin i^y sin 0y sin <£y + cos i//y cos <£y^ 

+ ^w + py 2 — qx^sin 0y cos 0y 

Wj = (u + qz x - ry x ^(cos \p z cos sin 0 Z + sin ^Z sin ^z) 

+ (v + rx y - pzy^sin \p z cos (p z sin 6 Z - cos \p z sin (/> z ) 

+ (w + py z - qx^cos cp z cos # z 
Angular velocity. - 


= p cos $ x cos i// x + q cos $ x sin i// x - r sin $ x 

q. = p ^cos sin 0 y sin 0 y - sin i// y cos 0 y ) 

+ q^sin 0y sin i^y sin 0y + cos i//y cos ^y)+ r s* 11 0y cos ^Y 

= p^cos \p z cos <p z sin 0 Z + sin \p z sin cp^j 

+ q^sin \p z cos (p z sin 6 Z - cos ^ z sin <£ z ) + r cos 6 Z cos cp z 

Transfer From Instrument Axes to Vehicle Reference Axes 


(1-7) 


(1-8) 


(1-9) 


If the axes along which data are measured are not orthogonally alined, the following 
procedure must be used to transfer instrument readings to the vehicle reference axes 

at the c.g.: 

Accelerations are given by 


gA X,cg 

gA Y,cg 

- [j'nonorthogonal] 

s A x,r 

gA Y,i 

+ 

K 2 

jS A Z,cg_ 


_ gA Z,i_ 


k 3 



(1-10) 
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SECTION I - Continued 


angular velocities, by 



P 

- fr l" 1 


N 



q 

[j nonorthogonalj 

^i 




r 


_ r i_ 



and linear velocities by 





u 


a i 


K 4 


V 

P nonorthogonalj 

Vi 

+ 

k 5 


w 


Wi 


K 6 


where 

K x - x x (r 2 + q 2 ) - y x (pq - r) - z x (rp + q) 
K 2 = -x y (pq + r) + y y (p 2 + r 2 ) - z y (qr - p) 
K 3 = -x z (pr - q) - y z (qr + p) + z z (q 2 + P 2 ) 
k 4 = ry x - qz x 
K 5 = pz y - rx y 
K 6 = qx z - py z 


( 1 - 11 ) 


( 1 - 12 ) 


(1-13) 


In equations (I- 10) to (1-12) [rnonQ^hogonalJ 1 is the inverse of a transforma- 
tion matrix made up of the sines and cosines of nine Euler angles, one set i// for 
each of the three instruments. This transformation matrix is defined as 


cos # x cos i// x 


cos 0 X sin i// x 


-sin # x 


f nonorthogonalj 


cos i// Y sin s i n 
- sin x^y cos <£ Y 


cos cos s i n e z 
+ sin i// z sin 


sin \1>y sin 0 Y sin sin 0 Y cos 0 y 

+ cos x^y cos 4>y 

sin xp^ cos 0 Z sin 6 z cos 6 z cos 0 Z 
- cos xp^ sin $ z 


(1-14) 


-1 


Except under certain simplifying assumptions, the inverted matrix |T nonor thogonal) 
needed to solve these equations is too long and involved to write out; however, the matrix 
inversion can be performed fairly easily on a digital computer. Equations for two cases 
in which the matrix can be inverted are given in appendix A as equations (A-5) to (A-16). 
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SECTION I - Continued 


CORRECTIONS TO ANGLES OF ATTACK AND SIDESLIP 


Angle of Attack 

For a vane -type sensor, the true or corrected angle of attack (referred to the vehi- 
cle reference axes at the c.g.) is given by 

a = - a a - q? u + a Y - Kq - + Aa (1-15) 

where Aa is the position error, the error due to the location of the sensor in the flow 
field of the body, and the other corrections are those due to boom bending, misalinement , 
and flow -field effects (see the list of symbols). 

Equation (1-15) is for an angle -of -attack vane mounted on an instrument boom extend- 
ing from the body or a wing tip but generally applies to any type of sensor. Other types 
of sensor and methods of calibration are discussed in reference 5. The importance of the 
different types of error is discussed in reference 3. 


The angle is the correction for vehicle rotation and is given by 

_ 2x ql 2y p£ 

a <i-TzvZ ' l 2 vZ 


(1-16) 


The other errors are determined by calibration as discussed in reference 5. 


Angle of Sideslip 

The equation for corrected sideslip angle is 

£ = + i3 a + 0 U - i3v + + /3 d + A/3 

where the correction for vehicle rotation is given by 

q _ 2x r£ 2z pi 

" £. 2Voo H 2V 00 


(1-17) 


(1-18) 


DETERMINATION OF FREE -STREAM MACH NUMBER 


For flight tests of high-speed aircraft or missiles, the following relationships are 
used to determine Mach number from onboard measurements of impact pressure q c 
and static pressure p: For subsonic conditions (M < 1.0), 


p 



JL. 




= (1 + 0.2M) 7 / 2 - 1 


(1-19) 


For supersonic conditions (M > 1.0), the equation is modified to include the loss in total 
pressure behind the shock wave and becomes 
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SECTION I - Continued 


<*c _ y + 1 



. 1 = !.2m 2 Lmsm! f 2 


- 1 


,5.6M^ - 0.8, 


(1-20) 


The indicated Mach number M is determined from tables or plots, based on equa- 
tions (1-19) and (1-20), of q c /p as functions of Mach number. 

DETERMINATION OF AERODYNAMIC FORCES AND MOMENTS FROM ' 

ACCELEROMETER AND RATE -GYRO READINGS 


The components of the total force acting on the flight vehicle can be determined by 
multiplying the corrected accelerometer reading (i.e., corrected for displacement and 
misalinement) by vehicle weight. In coefficient form 



WA 




X.cg 


q S 


WA 


q S 


WA 




q S 

OO 


( 1 - 21 ) 


The total force thus determined includes the static and the dynamic aerodynamic forces, 
the engine thrust, the jet damping force, and the reaction and other control forces. It 
does not include the components of vehicle weight, however, since gravitational effects 
appear as accelerations that are measured directly by the accelerometer; that is, the 
products of the weight and the corrected accelerometer readings are 

WA X,cg = W(a + wq - vr + g sin 6 g ) = ^ F x | 


WA Y, cg = W (y + ur - wp - g cos sin = ^ F Y / (1-22) 

WA Z, cg = W (^ + vp - uq - g cos 0 g cos 0 g ) = ^ F Z 
where the sine and cosine terms are components of vehicle weight. 


The total moments are determined from the rate -gyro readings according to the 
equations 
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SECTION I - Continued 


^M x = pl x - qlxY ' rl xz + q(rl z - pl xz - qIy Z ) ‘ r ( qI Y " rI YZ " pI Xy) 

= qly - fryz “ P*XY + r ( pI X " qI XY ~ rI Xz) " p ( rI Z " pI XZ “ qI Yz) ^ 

^M z = rl z - Plxz " qI YZ + p ( qI Y “ rI YZ " pI Xy) " q ( pI X " qI XY " rI Xz) 

The moment coefficients are 



(1-23) 

4 


(1-24) 


(Although b and c are commonly used to nondimensionalize the moments, these lengths 
are arbitrary and any convenient length can be used.) 

Methods by which the total forces and moments are broken down into trim, static, 
and cross -coupling components so that longitudinal and lateral coefficients can be deter- 
mined are quite complicated and beyond the scope of this paper. The simplified method 
for determining basic static- and dynamic -stability derivatives from oscillatory flight 
motions is developed in appendix D. In this and in most methods, linear systems are used 
to represent transient flight data; these linear systems are obtained by assuming small 
angles, constant aerodynamic coefficients, constant free -stream conditions, rigid-body 
mass and inertial characteristics, and separation of the longitudinal and lateral modes. 

For data that do not fit these limitations, exact methods of simulation based on the equa- 
tions of motion for six degrees of freedom are used. Other methods that are now becom- 
ing popular involve parameter identification in which linear and quasilinear estimation 
techniques are used. (See, for example, refs. 2 and 6.) 


DIRECT FLIGHT MEASUREMENTS OF LIFT AND DRAG 


Lift and drag components for an airplane can be measured directly by using flight - 
path accelerometers which differ from body -fixed accelerometers in that they are oper- 
ated by a vane or pressure sensor that rotates the sensor into or normal to the stream 
direction. Accelerometers measuring parallel and normal to the flight path measure 
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vehicle thrust minus drag and lift, respectively, if the component of thrust in the lift 
direction can be assumed negligible. In equation form the lift and drag are 


t d - D = -wa d 

L = WA l 


(1-25) 


where Tp is the component of thrust along the flight path. The accelerations should be 
corrected for the usual boom -type errors as described in correction with equation (1-15) 
for the angle of attack measured from a boom. Details of this flight -testing technique 
are given in reference 7. 
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SECTION II 


TRANSFER OF AERODYNAMIC FORCE AND MOMENT COEFFICIENTS 
AND DERIVATIVES TO ANOTHER REFERENCE CENTER 

4 

The equations in this section are the general forms for transferring aerodynamic 
force and moment coefficients and stability derivatives from a coordinate axes system 
with origin at the c.g. to a parallel axes system with origin at a point x,y,z away from 
the vehicle c.g. (Simplified forms obtained by assuming zero angles of attack and side- 
slip and neglecting aerodynamic cross derivatives are given as equations (A-17) to (A-43) 
in appendix A.) In these general equations, unprimed coefficients are referred to the c.g.; 
primed coefficients are referred to the point at x,y,z shown in figure 2. 

The transformations are given for both systems of variables used in aerodynamics. 
Equations for the Q!,j3,V,p,q,r system are given by equations (II - 1) to (11-40); those for 
the u,v,w,p,q,r system are given by equations (11-41) to (n-76). (The relationships 
between the two sets of variables are given in table II.) The equations are derived in 
appendix B. They are general in that no assumptions are made as to angle of attack or 
sideslip and in that all transfer distances, lateral as well as vertical, are included; how- 
ever, they are still not complete in that, for the purposes of these transformations, 
second-order derivatives with respect to time are assumed negligible and are omitted. 
Also, the transformations given for the static forces and moments (derivatives with 
respect to a or 0) are the simplified forms that apply only when the body is not under- 
going any significant rotation and p, q, and r are essentially zero. The transforma- 
tions for static -stability derivatives that apply when there is significant vehicle rotation 
can be derived as indicated in appendix B. 

TRANSFORMATIONS FOR a,p,V DERIVATIVES 
X-Axis Force Coefficients and Derivatives 


,T - c 
X “ u x 

(H-l) 

l T — p 

Xq, - L x a 

Cn-2) 

% = % 

(II -3) 

X A = C Xq, 

(II -4) 

x • = c x • 

X /3 p 

(H-5) 
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SECTION II - Continued 


Cy =Cy - 2Cv 2. cos - a + 2Cv_ (f COS j3 + J sin a sin 0 
Xn a d £ cos )3 A /3\£ 


£ 


+ 4Cx( j sin |3 - j sin a cos /3 


(II -6) 


C^ = Cv + — - — Cy cos a + | sin a)+ 2C X/3 (f cos a« - y sin alsin j3 

Xq ^ cos|3 A «\! £ j A 0\£ £ / 


+ 4Cx(j sin a - ~ cos a J cos /3 

iX- x. 


(II -7) 


y sin a 


c’ =Cy - 2C V \ — — - 2C X J £ cos j3 + Z cos a sin /3 
X r A r £ cos p A j3\£ £ 


+ 4C x (X cos ce cos /3 - j sin /3 


(II -8) 


Y-Axis Force Coefficients and Derivatives 


Cy - Cy 
C Ya " °Y a 




% = C Y 
Cy . = Cy 

X OL 


Cy . = Cy . 


(II -9) 


c' = Cy - 2 Cy I £OS_a + 2Cy (f cos j3 + f sin a sin 0 
Y p *p * ol £ cos 0 x j3\£ £ 


+ 4Cy(^ sin /3 - \ sin a cos 0 


(II -10) 


2C- 


Cl - Cy + — [ £ cos a + % sin a ] + 2C V ( J cos a - £ sin aisin /3 

Y q *q cos /3U 1 / Y j3\4 X / 


Q? / x z 

1 £ COS Q! + 


X 


+ 4Cyf j sin a - j cos O') cos j8 


(II -11) 
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- <- - 2C,r £ si — - 2Cv (t COS 0 + \ cos at sin 0 

Y r Y r Y ot £ cos 0 Y 0\£ * 

+ 4Cy(| cos a cos i 3 " | sin i 3 ) 


(n-12) 

* 


Z-Axis Force Coefficients and Derivatives 


C Z = C Z 




c' = c z 

"OS 05 


% = c zp > 

cU =c z . 


■a 


0! 


C’ 7 . = c z . 

Z /3 Z |3 J 


C’„ -C y _-2C K \ + 2c Zp[j cos (3 + j sin a sin 0 


z p £ cos ^ 


(11-13) 


+ 4C Z ( j sin 0 - sin at cos 0 


(H-14) 


2C 


C’ = Cry + 


Z 


z q z q cos 0\i 


cos at +| sin atj + 2C Z ^| cos at - | sin atjsin 0 


+ 4C-z(“ sin at - j cos atj cos 0 


(11-15) 


c* = Cry - 2 C ry Y - 2C 7 £ cos 0 + £ cos at sin 0 

c Z r ^ z r z at e cos 0 Z 0U * 


x . 


+ 4C Z [ j cos at cos 0 - - sin 0 


(11-16) 


X-Axis Moment (Roll) Coefficients and Derivatives 


= C Z + ! C Y-7 C 


l ~Z 


C i a = C la + \ Cy <* ' t C Z« 


(11-17) 

(11-18) 
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V'Vi CY «-f% 

°lp = % + I Cy /3 " f Cz |3 


% = <V I C Yp ' f % ' ¥ ^ (% + I C Yo; - f C Z a 


+ 2\j cos /3 + | sin a sin ^ + | C Y/3 - | C Z/3 


+ 4(| sin ft - | sin a cos ftj\C t + J Cy - J C z 


2 /x 


C Zq = % + f % - £ C Zq + 35^ cos “ + f sin a )[ C l„ + f C Y, 


(x & 1 & 


+ 2 sin /3(9 cos a - - sin a| (c> + - Cv„ - ^ Cry 
\£ J f l\ l ft i Y ft £ ^Zp 


+ 4 cos 0(| sin a - | cos ajfCj + 1 C Y - ^ C z 


C i r = %+! C Y r ' f C Zr - 2 | + f C Yq , - | C Zff 


X 


- 2 lf cos /3 + | cos a sin ft) (c lp + £ C Y(3 - f C Zj3 


+ 4(| cos a cos ft - | sin ft) (Cj + | C Y - £‘C Z 


Y-Axis Moment (Pitch) Coefficients and Derivatives 


C m " C m + f C Z "f C X 


C* =C + — P - — P 
m Q! m CX £ Zq; £ Xq; 


(H-19) 

( 11 - 20 ) 

(11-21) 


(11-22) 



(11-23) 


(11-24) 


(11-25) 


(11-26) 



SECTION n - Continued 


C m /3 = C ™0 + I Cz j3 " I ° X (3 
c m^ = C m^ + J C Z A " J C X^ 

C m^ = + J C Z^ - J C X^ 


C m D = C mp + f C Zp "f C X p - 2 fl|f ( C m « + fC Zo! "CXq, 


+ 2( j cos 0 + 1- sin a sin P)[C m p + J c Zp " J C X^ 


+ 4 | sin 0 - ?■ sin a cos 0 lC m + jC z - rC x 


C f - Crv, +-C7 --C Y + — 5-rff cos a + | sin a)(c m + y C z 
L mq Ly niq £ z q £ Xq cos 0\£ C. ) \ a l a 


+ 2 sin 01 j cos o? - j sin aj (c m £ + j C z ^ - j C x ^ 


+ 4 cos /3^ sin a - | cos aj [c m + | C z - | C x 


c m r “ C m r + J c Z r “ ? C X r " 2 f cos 3 V Cm « + 8. Cz « ~ 8 ° x oi 


L r ft £ A r £ cos 0 


- 2f| cos 0 + | cos a sin 0 J \C m ^ + j C z ^ - | C X/3 


+ 4(£ cos a cos 0 - x sin 0] lc m + x C Z ~ T C X 


x 


Z-Axis Moment (Yaw) Coefficients and Derivatives 


C'n = C n + I C X-I C Y 


c n a = C nce +f C Xo! -f CY 


f 


(11-27) 

(U-28) 

(11-29) 

(11-30) 

- ! c x a ) 

(11-31) 

(11-32) 

(11-33) 

(11-34) 
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SECTION II - Continued 


C Y = Cv 

A w A w 

(11-43) 

C Xp = C Xp + l C X V - \ c x w 

(11-44) 

C Xq = C Xq - f C X U + f C X W 

(11-45) 

C X r = c x r + 1 c x u - f c x v 

(11-46) 

Y-Axis Force Coefficients and Derivatives 

C^r = Cy 
Y u Y u 

(11-47) 

cV - c v 

x v *v 

(11-48) 

C Y W = °Y W 

(11-49) 

C Y p = % + f C Y V ' I C Y W 

(11-50) 

C V q - C Yq ' I C Yu - f C Yw 

(11-51) 

C v — Cy + r Cy - t Cy 
Y r Y r i *u & x v 

(n-52) 

Z-Axis Force Coefficients and Derivatives 

C’ = C v 

(11-53) 

cV = c z 

Zy ^V 

(11-54) 

C 17 = C 7 

(11-55) 

C Zp = C Zp + f C Z V c z w 

(E-56) 

c ’z q = c z q -I c z u + I c z w 

(11-57) 

C Z r = C z r + | C Z U ■ I c z v 

(11-58) 
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X-Axis Moment (Roll) Coefficients and Derivatives 

C 'lu = % + f C Y U -f C z u 

c ; v -% + l c y v - ! c z v 


C ^w = % + l C Yw-f C Z w 


C L = C ln + T C Yn - f C Z + + f C Yv C Z 


p £ Y p l i l 4 V £ X V £ "v 


- T C Z +f C Y -T C Z 
£ \ a Y w £ 


c: = 


= C 7 +;c Y -£C 7 -tC, +JC Y -£C 7 
£q ^q l Y q fL /<q £ \ i Y u 1 


x 


+ !I c *w + ! c yw-! c zw 


c! = C 7 + 7 CV -^-C-z 4 c, + - C v 4 c 7 

£r £ Yr £ z r £ Wu 4 Y u £ z u 


-9 C; + x C v -*■ C 7 
M W £ Y v £ z v 


Y-Axis Moment (Pitch) Coefficients and Derivatives 

C m u = S + | C Zu - f C Xu 


c m v = c mv + f C Zy - § C Xy 


Cm,.,= Cm w +f C Zw -f C x 


w 


w 


X 


'm t 


- C mp + J C Zp - | C Xp + | ^C my + j C Zy - j C x 


- n c m w + f C Z W - I C X W 


C m q = Cm q + I C Zq ‘ f % - f K - f ° Zu - J C Xu 


+ fl Cm w + f c z w " I c x w 


( 11 - 59 ) 

. (n-60) 

( 11 - 61 ) 

\ 

( 11 - 62 ) 

( 11 - 63 ) 

( 11 - 64 ) 

( 11 - 65 ) 

( 11 - 66 ) 

( 11 - 67 ) 

(11-68) 

( 11 - 69 ) 
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c m r - c m r + f C Z r " f C X r + f ( C m u + f C Z U “ f C X U 


- ~(C m + — Cy --Cv 

0 \ niy £ Zy £ Ay 


Z-Axis Moment (Yaw) Coefficients and Derivatives 


Ck u = Cn u+ f c Xu-! C Y u 


Cn = C^ C x ~ ~r C Y 

n v n y £ Ay £ ly 


Cn = c n + I C x ~yCy 
n w w i A w a Y 


C n p = % + 7 C Xp -I C Yp + l( C n v + f C x v " J Cy v 


- 1 l c n w + J C X W - J C Y W 


C nu+ |c Xu -fc Yu 


+ ? (C n __ + 1 Cv -?CS 


w ’ £ ~" x w £ *w 


Cn = Cn + T C X ~ "T Cv + “r [Cn + ~ C x - — C x 
n r n r £ £ x r d \ n u £ A u £ Y u 


- - C n +^Cv - - Cy 
£ \ n v # A v a *v 


( 11 - 70 ) 

£ 

( 11 - 71 ) 

( 11 - 72 ) 

( 11 - 73 ) 

( 11 - 74 ) 

( 11 - 75 ) 

( 11 - 76 ) 
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SECTION III 


TRANSFER OF AERODYNAMIC FORCE AND MOMENT COEFFICIENTS 
AND DERIVATIVES TO A ROTATED AXES SYSTEM 

The equations in this section are the general forms of the axes transformations for 
aerodynamic force and moment coefficients and stability derivatives that can be used to 
transfer between any two of the five axes systems used in aerodynamic analysis. (Par- , 
ticular forms for transfer between body and wind-tunnel stability axes are given in 
appendix A.) Transformations are presented for both a, j3,V derivatives and u,v,w 
derivatives. (For the relationships between the two sets of derivatives, see table II.) 

The equations are developed in terms of a generalized angle of attack A and generalized 
sideslip angle B; transformations between any two axes systems can be obtained by sub- 
stituting specific angles for A and B in these equations. (See tables III and IV.) 

DESCRIPTIONS OF AXES SYSTEMS 

The five axes systems considered are those shown in figure 3. While three of 
these systems, body, principal, and wind axes, are clearly defined in the literature, there 
is some confusion concerning the definition of the stability axes. The stability axes 
described in some reports are vehicle or flight stability axes about which the equations 
of motion are written; in other reports they are wind-tunnel stability axes about which 
aerodynamic data are measured in the wind tunnel. The differences between the two are 
pointed out in these brief descriptions of the axes systems. 

Body Axes 

The orthogonal body -axes system is fixed within the vehicle with the X-axis along 
the longitudinal center line of the body, the Y-axis normal to the plane of symmetry, and 
the Z-axis in the plane of symmetry. This is the axes system about which aircraft 
instruments are usually mounted. Its main advantage in motion calculations is that vehi- 
cle moments of inertia about the axes are constant, so that the I terms can be omitted 
from the equations of motion. It is the logical system to which to refer velocities, accel- 
erations, and stability and control parameters in the study of aircraft handling qualities 
because the pilot T s orientation with respect to this frame is fixed. 

Principal Axes 

The principal axes are an orthogonal body -fixed system for which the products of 
inertia are zero. The X and Z principal axes lie in the plane of symmetry; the angle 
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SECTION III - Continued 

between the X body axis and the X principal axis is usually small so that in many 
cases the body axes can be assumed to coincide with the principal axes. 

Flight Stability Axes 

The flight stability axes (sometimes referred to as vehicle stability axes) are an 
orthogonal body-axes system fixed to the vehicle, the X-axis of which is alined with the 
relative wind vector when the vehicle is in a steady -state trim condition but then rotates 
with the vehicle after a disturbance as the vehicle changes angle of attack. This system 
is preferred in many stability studies because, as with other body -fixed axes, the moments 
of inertia about the axes remain constant and also because the motions defined are pri- 
marily those about the flight path rather than about body reference lines. 

Wind-Tunnel Stability Axes 

The wind-tunnel stability axes are the system about which most wind-tunnel data 
are obtained. For this system the X-axis is in the same horizontal plane as the relative 
wind at all times (fig. 3). The angle a between the X-axis of this system and the X 
body axes is variable. (It is a constant otQ for the flight stability axes.) This means 
that vehicle moments of inertia about the X-axis change. It also means that additional 
terms are required in the transformation equations for static -stability derivatives and 
for u,v,w derivatives when data are transferred to or from the wind axes or the wind- 
tunnel stability axes. These additional terms are designated (1), as in equation (III-7), 
for example. 


Wind Axes 

The wind axes are the system generally used in calculating motions of the vehicle 
as a point mass. The X-axis for this system is alined with the relative wind at all times 
so that vehicle moments of inertia about this axis change. As with the wind-tunnel sta- 
bility axes, additional terms, designated (1), are required in the transformations to or 
from the wind axes and either the body, principal, or flight stability axes, since the angle 
A between the X wind axis and the X-axis of either of these systems is variable. Also, 
since the lateral angle B between the X-axes is variable, there are additional terms, 
designated (2), as in equation (III-13), required in the transformations for some of the 
lateral derivatives between the wind axes and either of the other axes systems. 

NOTES ON USE OF TRANSFORMATION EQUATIONS 

In the transformations that follow, symbols such as Cx^Y’^Z and Cx’^Y’^Z 
are used in a general sense to designate coefficients and derivatives about corresponding 
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X r ?Y T ,Z f and X,Y,Z axes systems; specific designations to use with each axes system 
are given in table V. Two types of transformation are given: "Direct" which include 
those defined in table III and "inverse" which include those defined in table IV. Trans- 
formations between any two axes systems are obtained by selecting the proper angles for 
A and B from table III or IV and then the proper coefficient designations from table V. 
For example, in transferring from flight stability axes to body axes (a direct transforma- 
tion, according to table in), angle A is replaced by and angle B equals zero; the 
transformation for C x then becomes, by using equation (III-l) and table V, 


C X = C X,s cos a 0 ~ C Z,s sin a 0 




As a second example, in determining the C t about the wind-tunnel stability axes 
from derivatives given about the body axes (an inverse transformation according to 
table IV), equation (III-82) and tables IV and V are used to obtain 

%,wt = C h COs2 01 + Cn r sin2 01 + ( C Z r + C n p ) sin a cos a 

In the transformations for static -stability derivatives and also in the transforma- 
tions for u,v,w derivatives, the terms designated (1) are included only in transferring 
from either wind or wind-tunnel stability axes to either body, principal, or flight stability 
axes and not in other transformations. For example, terms (1) are included when the 
body-axes derivative C Xq; (eq. (III-7)) is determined from data given about wind-tunnel 
stability axes as in 


C x a - C Xa wt cos or - C Vwt Sin a - C Xwt sin a - C Zwt cos a 

but are not included when C x ^ is determined from data about flight stability axes as in 
C x a = C X a>s cos “0 - C Za>s Sin ckq 

The terms designated (2) are included only in transferring to or from wind axes. 

For example, in transferring from wind axes, the body -axes derivative (eq. (III-16)) 

is given by ^ 


C Ir = C ln cos a cos 0 - C m 
P 6 /3,w m /3,w 


cos a sin - C- sin a 
n j3,w 


" w cos a sin & ~ ^m,w cos a cos P 
but, in transferring from wind-tunnel stability axes, becomes 

= c 7 ^ cos a - Cj, ^ sin a 
^/3,wt n /3,wt 


If the need arises, derivatives with respect to Q!,j 8,V can be converted to deriva- 
tives with respect to u,v,w by using the relationships given in table II. 
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DIRECT TRANSFORMATIONS 


Static Force and Moment Coefficients * (Direct, Table in) 


c x 

- c x 

cos A cos B 

i 

O 

cos A sin B - C 

z sin A 

(m-i) 

§ 

Cy 

ii 

O 

X 

sin B + Cy 

cos B 



(III-2) 

C? Z 

- c x 

sin A cos B 

- Cy 

sin A sin B + C z cos A 

(m-3) 

o 

= C^ cos A cos B 

- C m cos A sin B - C n 

sin A 

(m -4) 

c’ 

^m 

ii 

O 

0«i 

sin B + C m 

cos B 



(m-5) 

n 1 

^n 

= C^ sin A cos B 

- C m sin A sin B + C n 

cos A 

(m -6) 


Static -Stability Derivatives (Direct, Table in) 

C X a = C x a cos A cos B - Cy^ cos A sin B - c Z a sin A 

- sin A cos B + C y sin A sin B - C^cosA (III -7) 

(l) 

where the terms designated (1) are included only in transferring from either wind or wind 
tunnel stability axes to either body, principal, or flight stability axes and not in other 
transformations . 

Cy = Cv sin B + C v cos B (III -8) 

Yqi A Q! ' ’ 

Clz = Cy„ sin A cos B - Cv sin A sin B + Cy cos A 
£ oi A ol 1 a 

+ C x cos A cos B - Cy cos A sin B - C z sin A, (III -9) 

(1) 


1 Equations for control derivatives and velocity derivatives have this same form; 
for example, 

C x = C x cos A cos B - Cy cos A sin B - Cy sin A 
^a °a °a c a 

^Xy ~ ^Xy cos A cos B - Cy-y cos A sin B - C Z y sin A 
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C? = Ci cos A cos B - C m „ ( cos A sin B - Cn sin A 

*• OL i (y Xil 0L OL 

- C t sin A cos B + C m sin A sin B - C n cos A 


Cm a = C^sin B + C mQ , cos B 


C^ “ c l a sin A cos B “ C m(X sin A sin B + C nQ! cos A 
+ Ci cos A cos B - C m cos A sin B - C n sin A 


Cl = Cv n cos A cos B - Cv n cos A sin B - C7 - sin A 
Xj3 */3 

- Px cos A sin B - Cy cos A cos B 


where the terms designated (2) are included only in transferring to or from wind axes. 

CL - Cy sin B + Cv Q cos B + Cv cos B - Cy sin B (III-: 

*/3 i3 P t 


Cry = C Y sin A cos B - C v sin A sin B + C 7 cos A 
z /3 A /3 1 /3 P 

- Cx sin A sin B - Cy sin A cos B 


C '1 = C ip cos A cos B - C m ^ cos A sin B - C n ^ sin A 
- Ci cos A sin B - C m cos A cos B 


= c lp sin B + Cmp cos B + cos B - C m sin B 
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C^ = C t sin A cos B - C m ^ sin A sin B + C n ^ cos A 


- Ci sin A sin B - C m sin A cos B 

j 

( 2 ) 


(III -18) 


Dynamic -Stability Derivatives (Direct, Table III) 


'X r 


(% 


n 


Cx cos B - fC Y + Cv_.Ysin B cos B + Cv sin^ B cos^ A 




+ c Z r sin2 A - (Cx r + C Zp ^ cos B - ^Cy r + C z ^' ; sin bJ sin A cos A (III-19) 




C Y cos 1 

Aq 


B - Cy p sin^ B + ^Xp - Cy^ s ^ n B cos B 


cos A 


(C 2q cos B + C Z p sin B^ sin A 


(III -20) 


C X r = ( C X r cos B - c Y r sin B ) cos2 A - (C Zp cos B - C Zq sin B) sin 2 A 


C Xp cos 2 B + C Yq sin 2 B - (Cxq + C Yp \ sin B cos B - C Z] 


sin A cos A (III-21) 


'Yp cos2 B " ^Xq sin2 B + (C X p - ^Yq) sin B cos B 


cos A 


- ^Cy r cos B + Cx r sin B^ sin A 


(IEE-22) 


c Y q = c Y q cos2 B + c X p sin2 B + (Cxq + Cy p ) sin B cos B 


(IB-23) 


"y r = ^ c Y r cos B + C Xr sin cos A 


^Yp c °s 2 B - Cxq sin^ B + ^C Xp - ^Yq) s ^ n B cos B 


sin A 


(HI-24) 


C’ Zp = (C Zp cos B - C Z q sin B^cos^ A - ^Cx r cos B - Cy r sin B^ sin^ A 

.2 D t C* r-< -i v-\ 2 


C X p cos * B + C Yq sin^ B - + Cy^sin B cos B - C Z] 


sin A cos A (HI-25) 
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C Zq = ( C Zq cos B + c Zp sin B^cos A 

+ c Xq cos2 B “ c Y p sin2 B + ^C Xp - Cy^sin B cos B sin A (III-26) 

c 'z r = C Z r cos2 A + c X p cos2 B + C Yq sin 2 B - + C Y ) sin B cos B sin 2 A 

+ [(C Xr + C Zp ) cos B - ^C Yj , + C Zq ) sin b 1 sin A cos A (III -27) 

“ 1 

C l p = C ^P C ° s2 B " { Cl q + Cm p) Sin B cos B + C mq sin 2 B cos 2 A + Cn r sin 2 A 

+ “( c Z r + c n p ) cos B + (C mr + Cj^sin Bj sin A cos A (III-28) 

^lq = ^iq cos2 B - C mp sin 2 B + ^C^ p - C mq j s ^ n B cos B cos A 

" ( C nq cos B + C np sin B) sin A (III -29) 

c Z r = (Cz r cos B - C mr sin B^cos 2 A - ^C np cos B - C^ sin B^sin 2 A 

+ C Zp cos 2 B + C mq sin 2 B - (C^ q + C mp ^sin B cos B - C nr sin A cos A (III-30) 

C m p = c m p cos 2 B - Cj q sin 2 B + (C ip - C mq ^ sin B cos B cos A 

~ ^Cin r cos B + sin B^sin A (III-31) 

C m q = C m q cos 2 B + C Zp sin 2 B + (C Zq + C mp ) sin B cos B (III-32) 

^m r “ ^C mr cos B + C^ sin B^cos A 

+ c m p cos2 B - C Zq sin 2 B + (C^ - C mq ^sin B cos B sin A (III-33) 
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^np (^np cos B - C^ sin B^cos 2 A - cos B - C mr sin B^sin 2 A 


C7 cos 2 B + Crri- sin 2 B - 


m q sin B - (C^ + C mp ^ sin B cos B - C n] 

Criq = cos B + C np sin B^cos A 

+ Qq cos 2 B - C mp sin 2 B + ^Cj p - C mq J sin B cos B 


sin A cos A 


sin A 


C n - C n cos 2 A + 


C lp cos2 B + c m q sin 2 B - (C lq + C mp ^ sin B cos bJ sin 2 A 


+ [( C ^r + Cn p) cos B " ( C m r + c nq^ sin Bj sin A cos A 


u,v,w Derivatives (Direct, Table III) 


X u “ |^X U cos2 B " ( c x v + C Y u ) sin B cos B + C Yy sin 2 B 
+ [j( C X w + C Z u ) cos B + ( C Y W + c Z v ) sin Bjsin A cos A 

\cos B ( C X S * n2 A COS B " C Y sin2 A sin B + C z sin A cos A^ 


cos 2 A + Cr, sin 2 A 


J w 


( 1 ) 


+ Cx cos 2 A sin 2 B + C v cos 2 A sin B cos B 


( 2 ) 


'X, 


cos A 


C X V cos B " C Yu sin 2 B + (C Xu - C y \ sin B cos B 

■' _ 

■ ( C Z V cos B + C Z U sin B ) sin A - ^(c x cos A sin B + C Y cos A cos b)cos 

( 2 ) 


] 


(III -34) 


(III -3 5) 


(III- 3 6) 


(m-37) 


J 

(III-38) 
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C X W = ( C X W cos B “ ^Y w sin B^cos 2 A - (C Zu cos B - C Zy sin B^sin 2 A 

+ [ C X U cos 2 B + Cy v sin 2 B - (^X w + Cyjsin B cos B - C z ^j sin A cos A 


cos A 

\ COS B 


(c x sin A cos B - Cy sin A sin B + C z cos a) 


( 1 ) 


+ sin A sin B (c x cos A sin B + Cy cos A cos B^ 


( 2 ) 


Cy u = [^Y U C0S ^ B " ^X V B + (C Xu " ^Yy) S ^ n B COS C0S A 


C y “ ^Yy cos 2 B + ^Xu ®i^ 2 B + ^Cx y + Cy ^ sin B cos B 

+ ^X cos B - Cy sin b) cos B^ 

( 2 ) 

C’ = (Cy cos B + Cv sin B^cos A 
Y w V 1 w A w J 


+ jcy u cos 2 B - C Xy sin 2 B + ^C Xu - C Y V ) sin B cos bJ sin A 
“ (^X cos B " Cy sin B^ sin A sin B 


( 2 ) 


(m-39) 


- ( C Y W cos B + C Xw sin B) sin A - ^C x cos B - Cy sin b) cos A sin B (III-40) 

( 2 ) 


(m-41) 


(IH-42) 


C Zu = ( C Z U cos B - C Zy sin B^cos 2 A - (Cx w cos B - Cy w sin B^sin 2 A 

+ |^x u cos 2 B + Cy v sin 2 B - (C Xv + Cy^sin B cos B - C z ^jsin A cos A 


- (C x cos A cos B - Cy cos A sin B - C z sin A^ 


sin A 
cos B; 


( 1 ) 


+ v (c x sin A sin B + Cy sin A cos b) cos A sin B 


(III -4 3) 


( 2 ) 
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c 'z v = ( c z v c ° s B + c z u sin B ) cos A 

+ fc Xv cos 2 B - C Yu sin 2 B + (C Xu - C Yy ) sin B cos Bj sin A 
- ^C x sin B + C Y c os b) sin A cos B , 


( 2 ) 

C' = C Zw cos 2 A + C Xu cos 2 B + C Yv sin 2 B - (C Xy + C Yu )sin B cos Bjsin 2 . 

+ [(C Xw + C Z U ) C0S B - (C Yw + C Zv )sin Bjsin A cos A 

cos A. 

+ ^Cx cos A cos B - C Y cos A sin B - C z sin ^ CQS ^ 

( 1 ) 

+ (^Cx sin A sin B + Cy sin A cos B^ sin A sin B 

_ _ ” ( 2 ) 


J 


c; = 


'U 


C^ u cos 2 B - + Cm u J s i n B cos B + C mv sin 2 Bj cos 2 A 

+ Cn w sin 2 A + [-(C, w + C nu )cos B + (C mw + C % )sin Bjsin A cos A 
+ (Cl sin A cos B - C m sin A sin B + C n cos A) 


sin A 
cos B 


(1) 

cy 

+ (C i sin B + C m cos B^cos A sin B 

( 2 ) 


C = 


C t cos 2 B - C mu sin 2 B + (C^ - C my ) sin B cos B 


cos A 


- (C ny cos B + C nu sin B)sin A - ^ sin B + C m co s B) 


cos A cos B 

j 


( 2 ) 


(III -44) 


(in-45) 


(III -46) 


(HI -47) 
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C Z W = (% C0S B ' C % Sin B ) cos2 A - (Cn u cos B - C % sin B)sin 2 A 

% eos 2 B + C mv sin 2 B - (C^ + C nlu j sin B cos B - C n J sin A cos A 


- sin A cos B - C m sin A sin B + C n cos A) 


cos A 


cos B 


7 


( 1 ) 


+ fa sin B + C m cos Bj sin A cos A sin B^ 


( 2 ) 


C m u “ \^ m u COs2 B " C lv sin2 B + ( C Z U " Cm v) sin B cos cos A 

- (Cm w cos B + C lw sin b) sin A - fccosB - C m sin b) cos A sin B 


( 2 ) 


C m v “ c m v cos2 B + Cj u sin 2 B + + C mu ^ sin B cos B 

+ ^ cos B - C m sin B) cos B 


( 2 ) 


C m w “ ( c m w cos B + sin B) cos A 


'm u cos* B - Cj 


l v sin2 B + (C^ - C m ^sin B cos Bjsin A 
- (c t cos B - C m sin B^ sin A sin B 


( 2 ) 


(HI -48') 


(III-49) 


(III -50) 


(III -51) 
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C n = ( Cn u C0S B " Cn v Sin B ) cos2 A “ (' c l w cos B " Cm w Sin B ) sin A 

+ [ Qu cos^ B + Cm v sin^ B - (pi v + Cm u ) s ^ n B cos B “ ^ n wj s ^ n A cos A 


/ _ . . _ _ . . \ sin A 

- (c i cos A cos B - C m cos A sin B - C n sin Aj cQ -“g 

CD 

+ {pi sin A sin B + C m sin A cos b) cos A sin B 

" ( 2 ) 


J 


c k v = ( c n v cos B + c n u sin B ) cos A 


*• [c iy cos 2 B - Cm u sin 2 B + (C iu - C mv )sin B cos bJ sin A 
- (Cl sin B + C m cos b) sin AcosB 


( 2 ) 


Ck w = C n w cos ^ A + 


q u cos 2 B + C my sin 2 B - (c ly + C mu jsin B cos B^sin 2 J 


(Cj w + C„ u )cos B - (C mw + C nv )sin B 


sin A cos A 


+ (Ci cos A cos B - C m cos A sin B - C n sin A) 

V L / cos B 

( 1 ) 


V. 


cos A 


+ (Ci sin A sin B + C m sin A cos b) sin A sin B ^ 

"( 2 ) 


(III -52) 


(III -5 3) 


(IB-54) 
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INVERSE TRANSFORMATIONS 


Static Force and Moment Coefficients 1 (Inverse, Table IV) 


Cx = ^C^- cos A + c’ z sin A^cos B + Cy sin B 
C Y " -(c' x cos A + C’ z sin A^sin B + Cy cos B 
Cy = -C^ sin A + C^ cos A 
Ci = (C^ cos A + C^ sin A) cos B + C^ sin B 
C m =s ~(c'i cos A + C^ sin A^sin B + C^ cos B 


C n = -C^ sin A + C^ cos A 


Static -Stability Derivatives (Inverse, Table IV) 

Cx„, = (Cy cos A + C^ sin A) cos B + CL sin B 
a \ A a J la 

- (^X sin A ” cos A ) cos B 


(1) 


Cy™ = -( C* cos A + CL sin A] sin B + C* cos B 
lot \ Xoi ) ia 

+ (Cv sin A - CL cos A) sin B 

j\. /Lx / » 


a) 


Cz,,, = -C’ sin A + Cl, cos A - C^ cos A - CL sin A 
‘-‘a aq, i x z 


' a 


J 


a) 


(m-55) 

■y 

(III -56) 
(III- 57) 
(III- 58) 
(III -5 9) 
(IH-60) 


(III -61) 


(HI-62) 
(III -6 3) 


* Equations for control derivatives and velocity derivatives have this same form; 
for example, 

Cy, - I CL cos A + Cy sin A )cos B + C v sin B 
6 a \ X 5 a 6 a 1 Y S a 

C Xy = (^Xy COS A + C Zy sin A V os B + C Yy Sin B 
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SECTION in - Continued 


C^ s ^Cj cos A + C' na sin A^cos B + C^ sin B 
- (cj sin A - Cn cos a) cos B 




( 1 ) 


Cm ffl = ~{ C l a cos A ^n^ ° iU ~ T 


A + C^ sin Aj sin B + C mof cos B 
+ (cj sin A - C^ cos a) sin B^ 


( 1 ) 


= - C J sin A + Cn a cos A - cj cos A - Cn sin A 


' n ot ' la 


( 1 ) 


(in-64) 


(III -6 5) 


(IH-66) 


cos A + cL sin A^cos B + C v ^ sin B 


Cv = ( Cv CUS Jr^ -t* V- ry O -111 r%. I ±J t v 

x j3 V x /3 z /3 / x /3 

" ( C X COS A + C Z sin A ) sin B + c y cos B 


( 2 ) 


Cv^ = -fc’ cos A + c'ry Sin A]sin B + Cy cos 

Y /3 V X /3 Z /3 ) 


B 


-£x 


cos A + c’ z sin A^ cos B - Cy sin B 


( 2 ) 

C 7o = -CL sin A + C™ cos A 
Z/3 X^ 

C Z/3 = fe cos A + sin A^cos B + C' mp sin B 


- (c^ cos A + Cn sin A^ sin B + c' m cos B 

( 2 ) 


J 


(m-67) 


(m-68) 


(m-69) 


(III -70) 
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SECTION III - Continued 


C m/3 - - (C^ cos A + C n/3 sin Aj sin B + C m ^ cos B 
- ^C^ cos A + C^ sin A^ cos B - C^ sin B 


(2) 


(in-7i) 


c n^ = - c ^ sin A + c 'np cos A 


(IB-72) 




Dynamic -Stability Derivatives (Inverse, Table IV) 

cos ^ B 


CV COS 2 A + C'ry Sin 2 A + ( CL 

Ap V Aj> jli 


+ CL \sin A cos A 
V) 


+ C 1 sin 2 B + 


q 


Cv + c’ V os A + fcl + C 1 \sin A 
Xq Y p J \ Y r Z q y 


sin B cos B (III -7 3) 


C Y = ( Cl cos A + Cl sin A^cos 2 B - ( Cl cos A + cl sin A\sin 2 B 


" TXq 


y 


■p 


Cl cos 2 A + Cl sin 2 A + ( cl + Cl ^ sin A cos A - C T 


'X 


P 


x r - -Zp 


sin B cos B (III -74) 


C X 1 


Cl cos 2 A - Cl sin 2 A - ^cl - cl ^ sin A cos A~| cos B 


'X, 


x p -z, 


J 


+ ( cl cos A - Cl sin A\ sin B 


(III -7 5) 

Cl cos A + Cl sin A^cos^ B - ( Cl cos A + Cl sin A^sin 2 B 
, Y P Y r / V Xq z q J 

- Cl cos 2 A + cl sin 2 A + ( cl + Cl \sin A cos A - Cl Isin B cos B (III-76) 
Xp Z r y x r Ap) iq 

cos 2 B + C^ cos 2 A + C^ sin 2 A + (c'x r + ^Zp^ s * n A cos aJ ® 


L ^ 

^Xq + ^Yp) CQS ^ + (^Y + ^Z ^) s * n Aj s i n B cos B 


(m-77) 


c Y r “ (^ C Y r cos A — Cy sin A^ cos B 


-Cl cos 2 A + cl sin 2 A + 
r 


ii i ry oui n T ( C V “ C ry J Si 

v Y -**p *- i rJ 


- Crv isin A cos A 


sin B 


(III -78) 
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SECTION ni — Continued 


C f 7 cos 2 A - c’ sin 2 A - (c' - CL \sin A cos A 

A r \ Ap z r y 


+ ( Cl 7 cos A - Cv sin A^ sin B 


cos B 


j q 




(m-79) 


: Zq = ^ C Zq cos A " C Xq sin A ^ cos B 


-C f z cos 2 A + Cv sin 2 A + (c\r - C 


P 


'X 


'X p - ~Z r 


sin A cos A 


sin B 


Cz r = C T Z cos 2 A + C^ sin 2 A - (ci 

I* V ^ 


X r + C Zp) sin ^ cos ^ 


'h ” I ~z 


ft 


C 7 _ = I c; cos 2 A + c 
p 

< T ™2 


P 

+ CL sin^ B + 


n r sin 2 A + ^C^ + cj^ sin A cos Aj cos 2 B 


C 7 _ + Cjjj cos A + (Cyy\__ + Cn J*) sin A 


'm 


r T ^ 


(III -80) 
(III -81) 


sin B cos B (III-82) 


c ln = ( C L cos A + C ^n sin A ^cos 2 B - (c' m ^ cos A + Cm r sin A) sin* B 


’ — [ v/i wo jr*. t un olil /i j 

^ V *-q u / 

C, cos 2 A + C^ sin 2 


7 


C lr = 


P 


c! cos 2 A - 


A + ( C ^r + C ^p) Sin A cos A “ c m q 


sin B cos B (III -8 3) 


C^p sin 2 A - ^C’^ - C^ sin A cos A 


cos B 


c m r cos A - c mp sin A^ sin B 


(III -84) 


C m „ = ( C’ cos A + C^ r sin a\ cos 2 B - 


V 


( c i q cos A + C Iiq sin A ) si 


sin 2 B 


C Z p cos ' Z A + C n r sin2 A + K r + CjJsin A cos A - C^ q 


sin B cos B (III - 8 5) 


^niq - Cmq cos 2 B + 


C’p cos 2 A + Cn r sin 2 A + ^ + cl^sin A cos A 


sin 2 B 


\ + C m p^cos A + ^C mr + Cnq^sin A 


sin B cos B 


(III - 86 ) 
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SECTION III — Continued 
Cm r = (Cm r cos A - C„ lp sin A^cos B 

cos 2 A + Cn p sin 2 A + - c[ lr J sin A cos A sin B 


c nr» - 


Cn cos 2 A - c! sin 2 A - 
i1m Iy 


( c ; P - c ».; 


sin A cos A 


cos B 


+ fcl cos A - Cj sin A) sin B 
Hq l q J 


= ( C ^q cos A " c * sin cos B 


+ |-Cn cos 2 A + c! sin 2 A 

1 n p l ~ 


q J 

n 


+ ^Cjp - Cj^ s ^ n A cos sin 


B 


Cn r = Cn r cos 2 A + cj p sin 2 A - (c[^ + C„ p ) sin A cos A 


u,v,w Derivatives (Inverse, Table IV) 

cos 2 A + C r-7 sin 2 A + ( cL + CL \sin A cos a] cos 2 B 

a w *uj 


Cx u = j^x u cos ^ A + s * n ^ A + (^"X w + C Z u ) sin A cos 

B + ^ c x v + C Y u ) cos A + ( C Y W + C Z V ) Sin 


+ Cv sin 2 
*v 


sin B cos B 


V ( C X S1 


z 


sin A - C^ cos A^ sin A 


J 


( 1 ) 


(c'^ cos A + C^ sin A^ sin B - Cy cos B 


cos A sin B 

j 


( 2 ) 


'X, 


= ( C X V 


cos A + C nr sin A^ cos 2 B - ^C^- cos A + C^ r sin A^ sin^ B 


u 


w 


.^2 


Cv cos 2 A + cly sin 2 A + ^C^ + Cly ^ sin A cos A - C^ 


L 


'X 


•u 


J w 


(ra-87) 


(hi - 8a) 


(m-89) 
(III -90) 


(m-9i) 


w J 


sin B cos B 


■v 


- ^C^- cos A + C ^ sin A^ sin B + Cy cos B 


cos B 
J 


(III -92) 


( 2 ) 
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SECTION III - Continued 


'Z u sin2 A - ( C 3 


( Cv - C ly ^ sin A cos A 
\ A u ^w/ 


= fCv cos 2 A - C 

k w r x w 

+ (Cy cos A - Cy sin A) sin B - (c^- sin A - C’ 5 
\ w u V v 

(1 


cos B 


+ [(C^ cos A + C’ z sin a) sin B - Cy cos bJ sin A sin B 


z cos A^cos A 

U) 


( 2 ) 


: ^Cy cos A + Cy w sin A^cos 2 B - ( c x v COS A + C Z V Si 


2 sin A^ sin 2 B 

:' x cos 2 A + C T Zw sin 2 A + ^C^ w + c’ z Jsin A cos A - C^ 


c; 


L V 


sin B cos I 


- ^C^ sin A - C z cos A) 


sin A sin B 


cos B 




( 1 ) 


[( C X COS A + C Z sin A ) cos B + c y sin B COS A Sin B 


Cv = C’ cos 2 B + 
x v Y v 


C’ cos 2 A + CL sin 2 A + fcL- + C ' 7 ) sin A cos A 
X u Z w V "u/ 


sur 


sin B cos B 


( C k v + C Y u ) COS A + ( C Y W + C Z v ) Sin A ] Si 

cos A + C’ z sin A^ cos B + Cy sin bJ cos B 


( 2 ) 


C v = ( Cl cos A - Cy sin A^j cos B 
Y w l Y w Y u J 


-Cl cos 2 A + Cl y sin 2 A + 


'X 


w 


J u 


( c x u - c ’z w , 


sin A cos A 


sin B 


+ ^C^ sin A - C z cos A) 


cos A sin B 


cos B 




( 1 ) 

( C X C0S A + C Z sin A ) cos B + C Y Sin ®] 

( 2 ) 


sin A sin B 


(III- 93 ) 


(III - 94 ) 


B 


(III -9 5 ) 


(III - 96 ) 
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SECTION III - Continued 


C z = 


CL cos 2 A - cL sin 2 A - ( CL - CL ^ sin A cos A 
/<u A w V A u 


cos B 


+ ^C’ z cos A - C^- sin A^ sin B + ^C^- cos A + C z sin a' 


sin A 
cos Bj 


( 1 ) 


C Zy = ^C z cos A - C v sin A) cos B 


'X 


■v 


~Cy cos 2 A + cLr sin 2 A + ( C^ - C^ l sin A cos A 
"" A w V A n / 


J u 


Z w Z w 


= C 7 .. cos^ A + cl sin 2 A 


sin B 


X 


u 


-(c; 


+ Cl, Vin A cos A 


X w T ^Z u/ 


- (cl cos A + sin A) cos ^ 
l x z ' cos B y 

oT 


c i = 

6 U 


c! cos^ A + cl sin 2 A + (c! + Cl ^ sin A cos A 

t u n w \ 1A uy 

+ C^ v sin 2 B + f(c; v + Cm u )cos A +(Cm w + Cn v )sin A~ 


cos 2 B 


sin B cos B 


+ (c£ sin A - C^ cos a') sin A 


J 


(1) 


C^ cos A + C^ sin A\ sin B - C^q cos B 


cos A sin B 
J 


(2) 


Ci y = (^1^ cos A + c n v sin A^cos 2 B - |cln u cos A + Cm w sin A^sin 2 B 


C l u cos2 A + c n w sin2 A + 


+ Cl ^sin A cos A - C* 


n u 


m. 


sin B cos B 


(c^ cos A + C^ sin Aj sin B - C^ cos B 


cos B 


.J 


( 2 ) 


(III-97) 


\ 

(III-98) 


(III -99) 


(III-100) 


(III-101) 
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SECTION m - Continued 


C *w = 


C m u - 


Cm v = 


C J W cos2 A ' C % Sl " 2 A ' ( C Z U ' C "w) Sin A COS A 


cos B 


- (^Cj nw cos A - C^ sin A^sin B - (c^ sin A - C^ cos A) cos A 

( 1 ) 


J 


(C* cos A + C^ sin A) sin B - C^ cos B 


sin A sin B 

J 


( 2 ) 


(Cm u cos A + CJh w sin A)cos 2 B - ^C^ cos A + Cn v sin A^ sin 2 B 


Cj cos 2 A + Cn w sin 2 A + f C^ + C^ sin A cos A - C^ 


“1 


(III -102) 


sin B cos B 


- (c^ sin A - C^ cos A^ 


sin A sin B 


cos B 




( 1 ) 


^ cos A + Cn sin A^ cos B + C^ sin B 


cos A sin B 


(III -10 3) 


_y 


( 2 ) 


CL cos 2 B + 


c * u cos 2 A + c„ w Sin 2 A + (c\^ + C^sin A cos A 
( C J V + Cm u )cos A + fc; nw + Cn v )sin aJ sin B cos B 
^ cos A + C 2 sin A^ cos B + C^ sin B 


sin 2 B 


cos B 


(III -104) 


( 2 ) 
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SECTION m - Concluded 


'm w - ( c m w ^ m u 


cos A - C' mu sin cos B 
+ - c l w cos 2 A + Cn u Sin 2 A + (cj u - Cn w )sin A cos A 


sin B 


( c i sii 


+ ( Cl sin A - Cn cos A^ cos ^ sin B 


o 


cos B 


a) 


(c^ cos A + sin aJ cos B + sin b| sin A sin B 


( 2 ) 


: n u = |^n u cos2 A - c i w sin2 A - ( c j u - C„ W J sin A cos aJ cos B 
+ cos A - c'i v sin AJ sin B + (cj cos A + sin aJ — 


sin A 


cos B 


y 


(i) 


Cn v = cos A - C\ y sin Aj cos B 

+ pCn u cos 2 A + Cj w sin 2 A + (c'^ - c' n ^ sin A cos A 

c n w = Cn cos 2 A + C ] sin 2 A - (c\ + Cl\ sin A cos A 

w w t-u \ tiff uy 


sin B 




cos A + C n sin A 


cos A 
cos B 


a) 


(in-105) 


(m-106) 


(m-i 07 ) 


(in - 108 ) 
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SECTION IV 


TRANSFORMATION EQUATIONS FOR MOMENTS OF INERTIA 


This section gives equations for transferring moments of inertia to a rotated axes 
system in a general form, similar to that used in section III, so that transformations can 
be made between any two of the five axes systems used in aerodynamic analysis. As in 
section ni, these transformations are defined as direct and inverse in tables III and IV, 
respectively, and can be performed between any two axes systems by selecting proper 
angles for A and B from these tables. 

These equations are derived from a general transformation of the form 



(IV-1) 


where X is the transformation matrix given in appendix A and is applied to the aerody- 
namic axes by setting = -B, 6 - A, and 0 = 0. 

The most commonly used forms of these transformations, those for transfers among 
body, principal, and flight stability axes, are given in appendix A. 


DIRECT TRANSFORMATIONS (TABLE III) 


I x = (lx cos ^ B + 2 Ixy sin B cos B + Iy sin^ b) cos 2 A + Ig sin^ A 



+ ( 2I XZ 

cos 

B - 

2Iyz sin b) sin A cos A 

(IV-2) 

J Y = 

Iy COS 2 

B + 

o 

I x sin B - 21-^y sin B cos B 

(IV -3) 

I Z = 

Ig COS 2 

A + 

(lx 

cos 2 B + Iy sin 2 B + 21-gy sin B cos b) sin 2 A 



- (2IXZ 

cos 

B - 

2 IyZ sin B) sin A cos A 

(IV-4) 

t t 

a xy 

“ i xy ( cos2 

B - 

sin 2 b) - ^I x - sin B cos B cos A 



- ( X YZ 

cos 

B -t 

■ ^xz ®i^ B^ sin A 

(IV -5) 
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SECTION IV - Concluded 


J XZ “ ( x xz cos B - Iy2 sin b)cos 2 A - ^I^z cos B - Iy Z sin B)sin 2 A 
“ (^X cos 2 B + Iy sin 2 B + 2Ixy sin B cos B - Ig) sin A cos A 

Iyz = (* yz cos B + Ixz sin B ) cos A 

+ [^Xy( cos ^ b ~ s * n ^ B ) " (lx - Iy) s * n B cos bJ sin A 


INVERSE TRANSFORMATIONS (TABLE IV) 


l X = ( r x cos2 A + I Z sin2 A " 2I XZ sin A cos A ) cos2 B + Iy sin 2 B 
“ 2 ^L^y cos A + Iyz sin B cos B 

Iy = Iy cos B + ^I^r cos 2 A + I z sin 2 A - 2I^z sin A cos a) sin 2 B 
+ 2 ^I^y cos A + Iyz sin a) sin B cos B 

I Z = *z cos2 A + sin2 A + 2Ixz sin A cos A 


I XY = 

\}XY 

cos A + Iy Z 

sin a) cos 2 B - 

■ ^I^y cos A ‘ 

f X YZ 

sin A) sin^ 


+ ( r x 

cos 2 A + I f z 

sin 2 A - 2I^- Z 

sin A cos A - 

Iy) sin B cos B 

l XZ = 

[*xz( 

cos 2 A - sin 

2 A ) + " l Z, 

) sin A cos aJ 

cos B 



+■ ^Iyz cos A - I-^y sin a) sin B 




!yZ = | 

ftz 

cos A - I^y 

sin a) cos B + 

r xz ( si “ 2 A - 

cos 2 

a) 


' ( : X 

- I T Z ) sin A cos A sin B 





B 


(IV-6) 
(IV -7) 

(IV-8) 

(IV -9) 
(IV-10) 

(IV-11) 

(IV-12) 

(IV-13) 
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SECTION V 


EQUATIONS OF MOTION FOR SIX DEGREES OF FREEDOM 

The equations presented here are the general forms that include the variables „ 
likely to be of interest in computing motions in the vicinity of the earth (moon and sun 
perturbations are ignored). The terms of each equation are grouped so that various 
effects (for example, oblateness of the earth) can be accounted for by adding or omitting 
certain terms. Linearized equations and the wind-axes equations for a point mass are 
given in appendix A. 

The general equations apply to any of the five systems of vehicle reference axes 
shown in figure 3. The Euler angles %,0g,0g are referred to the gravity -axes sys- 
tem, with origin located at the surface of the earth, which rotates with the earth as shown 
in figure 4. (The relationship of these gravity axes to the vehicle reference axes is 
shown in fig. 5.) 


FORCE EQUATIONS 

Equations for the forces along the X,Y,Z axes are given in general form as equa- 
tions (V-l) to (V-3), respectively. Equations for specialized cases can be obtained from 
the general forms as follows: 

(1) For an oblate earth, equations (V-43) to (V-45) of auxiliary equations are used 
in place of term (2) in X,Y,Z force equations, respectively. 

(2) If mass of vehicle is constant (zero thrust is also implied), terms (3) to (5) are 
omitted. 

(3) For flight outside the atmosphere, terms (6) to (8) are omitted. 

X-Axis Forces 


m(u + qw - rv) + mg sin 0 g = T x - m(q d z “ r d y), + . F r,X. 
's ^ J V ^ S' v — ^ '-s 

(1) (2) (3) (4) (5) 


+ s( c x>0 + c X q! « 4- c x v ^ 

( 6 ) 


+ fop? + Cx p Jfc + °Xr WZ; + , Cx 6a 6a + j + ^ 


(7) 


( 8 ) 


(9) 
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SECTION V - Continued 


where the different terms are defined as follows: 


(1) Mass times acceleration 

(2) Component of vehicle weight (for a more nearly complete weight component, 

one that includes earth -oblateness effects, see auxiliary eqs. (V-43) to (V-45)) 

(3) Primary rocket thrust (see auxiliary eqs. (V-36) to (V-39)) 

(4) Jet damping force 

(5) Reaction control force * 

(6) Basic aerodynamic forces 

(7) Aerodynamic cross -coupling terms 

(8) Aerodynamic control forces 

(9) Higher order terms 

Expansions of the aerodynamic forces and moments (terms (6) to (9) in force equa- 
tions and terms (11) to (14) in moment equations) neglect all aerodynamic partial deriva- 
tives with respect to rates of change of velocities and angles except those with respect to 
oi and /3. The forces are expanded in terms of the independent variables of,/3,V,p,q,r 
but could as easily be expanded in terms of the variables u,v,w,p,q,r. The X-axis aero- 
dynamic force, for example, could also be written as 




X - g pV °° S ( C X,o + C X U U + C X V V + C X W W + C Xp 2 y" + C Xq 


_aL 

2Voo 


Cx r + ^ c Xs 6 + Hi S her order terms^ 


Y-Axis Forces 

m(v+ ru-pw^ - mg cos 6> g sin = Ty - ^m(r d x - p d z ^ + 


( 1 ) 


( 2 ) 


(3) 


(4) (5) 

+ |pv£s(c Yj0 + Cy^ + Cy y + Cy£ M- + Cy p + Cy r $ 


oo, 


( 6 ) 


+ C Y a a +C ?aWZ + Cy q $k, + + ° 


(V-2) 


(7) 


(8) (9). 


where the terms (1) to (9) are defined as in equation (V-l). 
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SECTION V - Continued 
Z-Axis Forces 


m(w + pv - qu) - mg cos 8g cos (pg - - m(p d Y - q dx) + F r,Z 

l ^ ) \ ^ -J 

(1) (2) (3) (4) (5) 


ipv^s (c z>0+ c Za <* + Cz v + cz* 4 - C Zq -ai- 


2V, 


( 6 ) 


+ c +c *?Mz + Cz p + CZr *h* + CZ6 * 6e J \° 


(7) 


( 8 ) 


JJ ^ 
(9) 
(V-3) 


where the terms (1) to (9) are defined as in equation (V-l). 

MOMENT EQUATIONS 


Equations of motion involving moments about the X,Y,Z axes are given in general 
form as equations (V-4) to (V-6), respectively. Equations for specialized cases can be 

obtained from the general forms as follows: 

(1) If vehicle mass is constant, terms (2) to (4) and (7) to (9) are omitted 

(2) If principal axes are used, terms (3) to (6) are omitted 

(3) If vehicle has a plane of symmetry (X-Z plane) but principal axes are not used, 

terms (3) and (5) are omitted 


X-Axis Moment (Roll) 



= m 


-p(d Y 2 + d z 2 ) + d x (q d Y + r d Z ] 


(7) 



+ fpV ( 


2 S£ 


■1,0 - V * °zy sk * % Mz * C; P zk + 

( 11 ) 




+ o 
(14) 


(V-4) 
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SECTION V — Continued 


where the different terms are defined as follows: 

(1) > (5), and (6) moments of inertia times angular accelerations 

(2) to (4) and (7) jet damping moments 

(8) Moments due to main rocket thrust 

(9) Moment due to reaction control 

(10) Moment due to gyroscopic action of engine rotating mass 

(see auxiliary eqs. (V-40) to (V-42)) 

(11) Basic aerodynamic moments * 

(12) Aerodynamic cross -coupling terms 

(13) Moments due to aerodynamic controls 

(14) Higher order terms 

Y-Axis Moment (Pitch) 



= m 


( d x 2 + d z 2 ) + d y( p d x + r d z) 


(?) 


+ ^(T X d z - T z d x ) +M r?Y +M 


( 8 ) 


( 9 ) 


rm 

( 10 ) 


+ |pvjs^c m) 


+ C m _Q! + C 


a 


m vr 

v on 


+ c 


m 


al 
a2Voc 


+ Cm, 


q £ 


q 2 v c 




( 11 ) 


+ C 


m/ 


+ C 


m 


0 2V* 


+ C 


m 


Pi 
P 2V 0 


+ C T 


rl 
lr 2V„ 


- + Cj 


+ C 


m 6. 




( 12 ) 


+ 0 


(13) 


(14) 


where terms (1) to (14) are defined as in equation (V-4). 


(V-5) 


1 See statement following equation (V-l) about terms (11) to (14) in expansion of 
aerodynamic forces and moments. 
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SECTION V - Continued 


rl z + pq( Iy - I X ) + ri z - qi Y Z, - pi X Z, + (<1 2 ‘ P 2 ) X XY ‘ <P r ^YZ, + , (qr ~ ^Z , 

(5) (6) 


( 1 ) 


(2) ' (3) (4) 


= m 


- r ( d x 2 + d Y 2 ) + d z(p d x + q d Y^ - ( T x d Y J t y d x) ; + + 


(7) 


( 8 ) 


(9) (10) 


S 


+ | P v£s£| C n;0 


+ c n ^/3 


+ c 


V , r . _g£_ . Cn _p£. 

n VvZ Cn £ 2Voo n P 2 V 0 
( 11 ) 


„ ai 
, a + Cn • Trrr - 
oi u a 2Voc 


+ C 


_ai 

”q 2V 0 


( 12 ) 


+ c ng 6 a + c ng 6 r ] + 0 

v a v .ZL-J) ^ 

(13) (14) 

where terms (1) to (14) are defined as in equation (V-4). 

AUXILIARY EQUATIONS 


(V-6) 


General equations that take into account various relationships and effects in develop- 
ing the equations of motion are given as equations (V-7) to (V-45). Equations for special- 
ized cases may be obtained as follows: 

(1) For flat earth, terms containing L and X are omitted 

(2) For nonrotating earth, £2 e , Vx> Vy, and Vg are omitted 

(3) For no surface winds, terms containing W x , Wy, and Wg are omitted 

Relationship Between Euler Angles and Angular Velocities 

The Euler angles specifying vehicle alinement with the gravity -axes system can be 
determined from the angular velocities p,q,r by the equations 

. r cos <f)„ q sin . , . \> . T 

\b„ i- + - + L cos ifscr tan Ogr + (^e + 'Vl sin L 

® cos 0g cos 0g ® 

+ cos L sin i//g tan 0g) (V-7) 


0g = q cos 


0g - r sin (j ) g - p sin ipg - (fi e + x) cos L cos 


(V-8) 
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SECTION V — Continued 
L cos i ft 


cos o 


g 


0g - P + Q. tan 0g sin <pg + r tan Og cos <pg + 

The inverse relationships are 

P = 0 g - i// g sin 0 g - L cos e g cos i p g - (n e + x)(cos L cos 0 g sin 0 g 
- sin L sin 0g) 


ff / • \ sin if/g 

& + \^e + Xjcos L - 


cos 0 , 


g 


(V-9) 


(V-10) 


q - e g cos (j ) g + $/ g sin 0 g cos 0 g - L (sin 0 g cos I0 g sin 0 g 

- (n e + x) cos L (sin e g sin if/ g sin <* g + cos 0 g cos 0 g ~) 
/_ • \ 


sin 0g - sin i// g cos <pgj 


- (ft e + x) sin L sin <pg 


r= 


~0g sin 0g + if/g 


cos 0g 

cos 0g cos 6g - L (sin 0g cos \f/g cos 0g + sin if/g sin 0gj 


9g - L (sin 0g 

- (^ e + x) cos L (sin Og sin \f/ g cos 0g - cos if/g sin 0g^ 

- (£2 e + x) sin L cos <pg cos 6g 

Vehicle Coordinates 

can be computed from the X,Y,Z axes velocity 
rating the equations 


Vehit 

Vehicle coordinates can be compu 
and vehicle Euler angles by integrating 

cos \pg cos Og + v(cos if / g sin Og sin <pg - sin if/g cos cpg'j 



Xg = U 


sin Og 

+ w(cos if/ g sin Og cos <pg + sin xf/g sin (V-13) 

y g = u sin i/'g cos Og + v (sin if/g sin Og sin 0g + cos i f/g cos 

+ w(sin if/g sin cos 0 g - cos if/g sin <£ g ) (V-14) 


Sg u sin Og + v cos Og sin <pg + w cos Og cos <fi 


J g 


(V— 15) 


r g = f? + y g 2 + z g 2 


(V-16) 
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Trajectory Parameters (See Fig. 6) 

Flight -path angle, longitude, and latitude can be determined from the equations 
-1 h 


y = tan' 


X = 


\j(x - R e £2 e cos L)^ + y* 


L g 


R e cos L 


- Qc 


• -y & 

L = —5. 

Re 

The range for a spherical earth is 
Range = R e cr 

where, for a total range of less than 50 000 feet, 

a = ^/(L - L 0 ) 2 + [(X - X Q ) cos L^) 2 
For a total range greater than 50 000 feet, 

-1 


<J = cos 
range for a flat earth is 
Range = \Jx, 


j^sin L 0 sin L + cos L 0 cos L cos 


(A - Xo)j 


2 + y 2 

g ^ y g 


(V-17) 


(V-18) 


(V-19) 


(V-20) 


(V-21) 


(V-22) 


(V-23) 


Angle of Attack, Sideslip, and Relative Velocity (See Fig. 7) 

Angle of attack, sideslip, and resultant relative velocity are related to components 
of velocity along the vehicle body axes by 

-1 


a = tan' 


u b 


-1 v b 

B - sin 1 — 

V 


v = % 2 + ^b 2 + w b 2 


u h = V cos a cos /3 = u. 


- W, 


‘b - — ^ — h - - b " v x,b ‘ vv x,b 

v b = V sin p = v b - V Y) b - W Y)b 


(V-24) 

(V-25) 

(V-26) 

(V-27) 

(V-28) 
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SECTION V - Continued 
w b = v sin a cos 0 = w fo - V z?b - W z?b 
where subscripts b; X,b; Y,b; and Z,b denote body-axes components. 


(V-29) 


Wind Corrections (See Fig. 8) 


Geostrophic (due to earth T s rotation).- The components of the geostrophic wind 


along the 

X,Y,Z vehicle reference axes are 



V x = R e ^e cos L cos % cos ^g 

(V-30) 


Vy = -R e ^ e cos L ( cos 0g sin ^g ” sin % sin ^g COS *%) 

(V-3l) 


Vz = R e £2 e cos L ( sin 0g sin % + sin % cos ^g cos *%) 

(V-32) 

Surface winds.- Components of the surface winds along the X,Y,Z vehicle refer- 

ence axes 

are 



W X = W (cos A f sin B ? cos Og cos i//g - cos A cos B cos 6g sin ^ g 



+ sin A t sin 0g) 

(V-33) 


Wy = w|cos A’ sin b' (sin 0g sin 9g cos \pg - sin i/'g cos •pg'j 



- cos A T cos B* (sin xpg sin 6g sin <pg + cos xf/g cos <figj 

- sin A t sin <fig cos 0gj 

(V-34) 


W z = wjcos A 1 sin B T (cos p g cos <fig sin 6g + sin \js g sin <£g^ 

- cos A ? cos B 1 (sin xpg cos 0g sin 6g - cos xpg sin 0g^ 

- sin A T cos <pg cos 0gl (V-35) 


Resolution of Engine Thrust and Torque Into Components 
Along Vehicle Reference Axes (See Fig. 9) 


Vehicle thrust and moment can be resolved into components along the vehicle refer- 


ence axes 


by the equations 

T = T 0 + — u + 
° 9u 

T X = T cos 0 rm 
Ty = T cos 


9T 

86 rpm 


6 rpm 


cos ^/ rm 
sin \p rm 


(V-36) 

(V-37) 

(V-38) 


64 



SECTION V - Continued 

Tz - -T sin 0 rm (V-39) 

• • • • 

and, if it is assumed that I rm = f 2 rm = 6 rm = i// rm = 0 , 

L rm = s i n ^rm " r cos ^rm s ^ n ^rm) (V-40) 

M rm = !rm^rm( r cos e rm cos \p Ym + p sin 6 > rm ) (V-41) 

N rm = Irm^rm(P cos 0 rm sin t/^ rrn - q cos 0 rrn cos ip Y (V-42) 


Components of Gravitational Acceleration Along X,Y,Z Vehicle Reference 

Axes With Earth-Oblateness Effects Included 

If the effect of earth oblateness on the acceleration due to gravity is considered 
important, terms (2) in the X,Y,Z axes force equations (V-l) to (V-3) should be replaced, 
respectively, by the following X,Y,Z axes weight components: 


GE 

mg x = m p-< 


J 2 (lf) ^ Sin L C ° S + f J3 (^) ^ Sin2 L ~ 1 ) cos L 


§^ 4 (lr) ^ sin ^ ^ _ ^ s * n l)cos l - . 


cos 6 g sin i//g 


1 - 5 sin 2 L - l) - 2 J 3 (5 sin^ L - 3 sin l) 

C* \t\ J \ R / 


8 


/a \ 4 / 

4 -£■) (35 sin 4 L - 30 sin^ L + 3j - . 
\R / 


sin 9 




(V-43) 


mg Y = m 


GE 

r2 


p 2 3 

J 2 ^ sin L cos k) + l'^3(^) (s sin^ L - l) cos L 


2 


f J 4 (lf) Sin3 L “ 3 sin L)cos L - . . . 


(sin if/g sin 6g sin <pg + cos i//g cos 0 g ^ 


1 - ir (3 sin 2 L - l) - 2 J 3 ^ (5 sin^ L - 3 sin l) 

_ & \ rt / \ R / 


a ^ 4 


- g- (35 sin 4 L - 30 sin 2 L + 3) - . 


sin $g cos 9g 


(V-44) 
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SECTION V - Concluded 



where 

GE = 398 601.2 ± 0.4 km3/sec 2 ((14 074 901.1 ± 14) x 10 9 ftVsec 2 ) 
a e = 6378.160 ± 0.005 km (20 925 721.8 ± 16 ft) 

and J 2 , J 3 , and J 4 are the second, third, and fourth zonal harmonics. The values 
from reference 8 are 

J 2 = 1082.7(1 ± 0.1) x 10" 6 

J 3 = -2.56(1 ± 0.1) X 10" 6 

J 4 = -1.58(1 ± 0.2) X 10’ 6 

In equations (V-43) to (V-45), which were derived from the gravitational potential 
given in reference 8, oblateness terms through the 4th order are considered and the 
earth’s longitudinal oblateness is neglected. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 16, 1972. 
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APPENDIX A 


SUMMARY OF FREQUENTLY USED FORMS OF AXES TRANSFORMATIONS 

AND EQUATIONS OF MOTION 

EULER ANGLE TRANSFORMATION BETWEEN 
TWO ORTHOGONAL AXES SYSTEMS 


For two orthogonal axes systems, an initial X,Y,Z reference system and an 
X' ,Y',Z' system obtained by rotating the initial system through Euler angles 0,0, 
(in that order), the transformations between the two systems are given by these 
relationships. 


Direct Transformation 

The transformation from the initial X,Y,Z system to the X ,Y ,Z system alined 
at Euler angles 0,0,0 with the initial system is 


(A-l) 


where 


x r 


X 

Y f 

= r 

Y 

_Z'_ 


Z_ 


r = 


cos 6 cos 0 

sin 0 sin 6 cos 0 
- sin 0 cos 0 

COS 0 cos 0 s ^ n 0 
+ sin 0 sin 0 


cos 6 sin 0 

sin 0 sin 6 sin 0 
+ cos 0 cos 0 

sin 0 cos 0 sin 6 
- cos 0 sin 0 


-sin 6 
sin 0 cos 6 

cos 0 cos 9 


(A -2) 


Inverse Transformation 

The transformation from the X ,Y f ,Z ! system back to the X,Y,Z system is 


r~ i 

X 


X T 

Y 

= r - 1 

Y r 

_Z_ 


X 


(A -3) 
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where 

cos 0 cos 0 sin 
+ sin 0 sin 0 

sin 0 cos 0 sin 
- cos 0 sin 0 

cos 0 cos 9 

TRANSFORMATIONS FOR ACCELEROMETER AND 
RATE -GYRO MEASUREMENTS 

These equations are simplified forms of the general inverse transformations for 
accelerations and angular velocities given in section I; however, they can be used in most 
practical applications to correct accelerometer and rate -gyro readings for displacement 
and misalinement. For cases that do not fit the assumptions made, the general forms in 
section I must be used. 

Case I 

Orthogonal Instrument Axes; No Restrictions on Misalinement Angles 

If X,Y,Z axes accelerometers or rate gyros making flight measurements are 
orthogonally alined then, even though the misalinements of the instruments with respect 
to the vehicle reference axes are large, their readings can be corrected by using these 
equations . 

Acceleration corrections . - 

= gA-^ cos 9 cos 0 + gAy ^(sin 0 sin 9 cos 0 - sin 0 cos 0) 

+ gAg i(cos 0 cos 0 sin 9 + sin 0 sin 0) + (r^ + q^)x^ 

- (pq - r)y x - (pr + q)z x (A-5) 

gAy, C g " g^ x ^ cos 9 sin 0 + gAy ^(sin 0 sin 9 sin 0 + cos 0 cos 0) 

+ gA z ^(sin 0 cos 0 sin 9 - cos 0 sin 0) - (pq + r)Xy 
+ (p 2 + r 2 )y Y - (qr - p)£ y (A -6) 
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gA Z eg “ " gA X i sin 6 + S A y i sin 0 cos 6 + % A z,i cos 0 cos e 

- (pr - q)x z - (qr + p)y z + (p 2 + q 2 )z z (A-7) 

Angular -velocity corrections . - * 

p = cos 0 cos 0 + q^cos 0 sin 0 sin 0 - sin 0 cos 0) 

+ r^(cos 0 cos 0 sin 0 + sin 0 sin 0) (A-8) 

q = p^ cos 0 sin 0 + q^(sin 0 sin 0 sin 0 + cos 0 cos 0) 

+ r^(sin 0 cos 0 sin 0 - cos 0 sin 0) (A-9) 

r = -p^ sin 0 + q^ sin 0 cos 0 + r^ cos 0 cos 0 (A-10) 

Case II 

Nonorthogonal Instrument Axes; Small Misalinement Angles 

If X,Y,Z axes accelerometers or rate gyros are not orthogonally alined but the 
misalinements with the vehicle reference axes are small, then these correction equations 
can be used. It should be noted that, in equations (A-ll) to (A-16), the angles 0,0,0 
are in radians. 

Acceleration corrections. - 


gA X,cg = g ( A X,i ' ^X A Y,i + e X A z,i) 

+ (q 2 + r 2 )x x - (pq - r)y x - (pr + q)z x 

(A-ll) 

gA Y,cg = g 0Y A X,i + A Y,i " l#> Y A Z,i > 

) - (pq + r)x y + (p 2 + r 2 )y y - (qr - p)z y 

(A-12) 

gA Z ,cg = g (" e Z A X,i + ^Z A Y,i + A Z,i 

) - (pr - q)x z - (qr + p)y z + (q 2 + p 2 )z z 

(A-13) 

Angular -velocity corrections.- 



P = Pi - ^ x q i + ®X r i 


(A-14) 

q = ^ Y p i + q i ‘ ^Y r i 


(A-15) 

r = -0 zPi + + r t 


(A-16) 
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SPECIAL FORMS OF TRANSFORMATIONS FOR AERODYNAMIC 
FORCE AND MOMENT COEFFICIENTS AND 
STABILITY DERIVATIVES 


Simplified Forms for Transferring Coefficients and Derivatives 
to Another Reference Center 

These transformations are the simplified forms obtained from the general equa- t 
tions of section II by assuming zero angles of attack and sideslip and neglecting aerody- 
namic cross derivatives. 


X-axis force coefficients and derivatives . - 
C X ~ C X 


Cx = Cy 

A C1! 

% = °Xa 

C Xq = C Xq + Y ° x a ' T ° x 


(A-17) 

(A-18) 

(A-19) 

(A-20) 


Y-axis force coefficients and derivatives . - 

Cl = & 


'Y = ^Y 

cV„ = c v 

Y /3 Y 0 


C Y . = Cy 

Y /3 

'Yp 


/3 


C V„ = C Yp + Y C Y ( 3 


= Cy-*Z C 


■0 


+ 


i£ 


c Y 


Z-axis force coefficients and derivatives.- 


c’z = c z 


•a 


= C, 


‘a 


C y . — C>7 . 


c% = + M C 

'q 


Zn - " z q 0!" 


4z 

T c z 


(A -21) 
(A -22) 
(A-23) 

(A-24) 

(A-25) 

(A-26) 

(A-27) 

(A-28) 

(A-29) 
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X-axis moment (roll) coefficients and derivatives. - 


c; = c 'i + |c Y -|c z 

(A-30) 

W !% 

(A -31) 

WN§ 

(A-32) 

% = % + I % + T Cl 0 + ^2 Cz a * °Y/3 

(A-33) 


(A -34) 

Y-axis moment (pitch) coefficients and derivatives.- 


C m = C m+fc z -|c x 

(A-35) 


(A-36) 

Cm^ = C m ^ + | C z ^ - | C Xdi 

(A-37) 

Cm q = Cm q - f C Zq - f CXq + 2 l C ma 2 * C Z(y 


2xz 0 4z ~ 4xz ^ , 4z 2 n 

— 9 - Cxq, " t Cm — 7 c z + - 7 - C X 

£ 2 £ £ 2 4 2 

(A -38) 

Z-axis moment (yaw) coefficients and derivatives.- 


Ck = c n+fc z -f C Y 

(A-39) 

% = Cn /3 ' I ° Y (3 

(A -40) 

Vv* 0 ** 

(A -41) 

Cip - % - f C Yp - ^ c ^« * f C "6 - f - % 

(A-42) 

ck r -c nr -fc Yr -fc n ^2£c Yii .fic„-Sc y 

(A-43) 
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Transfer of Coefficients and Derivatives From Body 
to Wind-Tunnel Stability Axes 

These equations are the most frequently used axes transformations. They convert 
the coefficients C x , C Y , Cg, etc., measured about body axes in the wind tunnel, into 
the coefficients about wind-tunnel stability axes C L , C Y , Cj-j, etc. They can be 
obtained from the general inverse equations (HI-55) to (IH-107) of section HI by replac- 
ing the angle A by ol and letting the angle B equal zero. 


Cj-j = -C x cos a - Cg sin a 


(A-44) 

C Y = C Y 


(A -4 5) 

Cl = C x sin a - Cz cos a 


(A -46) 

c Z,wt = c l cos a + c n sin ol 


(A -47) 

Cm = C m 


(A-48) 

^n,wt ” ~^l s ^ n a + ^n c °s ol 


(A -49) 

C|a = -C x cos a - C 7 sin a + C T 
u a> A « &ci L 


(A-50) 

Ct^) . = “C x . cos ol - 0 7 , sin a 
u oi a 


(A-51) 

% = % 


(A -52) 

Cy /3 = % 


(A-53) 

C L a = C Xq, sin « - c Zq, cos « * Cq 


(A -54) 

C L(J = C X d sin 01 ' C Z & cos a 


(A-55) 

%,wt = % cos “ + C n /3 sin a 

• 

(A-56) 

%,w t = % C ° S a + Cn $ Sin “ 


(A-57) 

%,wt = % cos2 a + c n r sin 2 a + (c lr + C„ p ) 

(Sin ol cos a 

(A -58) 

C ^r,wt = C h cos2 a - °n p sin 2 a + (c nr - C Zp ) 

1 sin a cos a 

(A -59) 

C « 6 ,wt = C h C0S “ + Cn 6 sin a 


(A-60) 
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c m a = c mo. 


(A-61) 

Cra d = C m<* 


(A-62) 

Cmq = Cmq 


(A-63) 

Cmg = 


(A-64) 

Cn 0,wt = sin “ + C n p cos oi 


(A-65) 

C n . = -C n . sin cx. + Cn • cos a 

(3,Wt p n j3 


(A-66) 

Cn p,wt = Cn P COs2 “ ' % sin2 “ + ( Cn r - %; 

! sin a cos a 

(A-67) 

Cn r,wt = Cn r cos2 “ + % sin2 “ - ( c l r + c n p j 

i sin a cos a 

(A-68) 

Cn 5,wt = ' C *6 sin 01 + C «6 cos “ 


(A-69) 


Transfer of Coefficients and Derivatives From 
Wind-Tunnel Stability to Body Axes 

These transformations convert coefficients about wind-tunnel stability axes Cl, 
C Y > etc., into coefficients about body axes Cx, Cy, C z , etc. They can be 
obtained from the general direct transformations in equations (III-l) to (III-36) of sec- 
tion III by replacing the angle A by a and letting the angle B equal zero. 


Cx = -Cjj cos ol + Cl sin a = -C^ 

(A-70) 

Cy = Cy 

(A -71) 

Cz = -C^ sin a - C L cos a = -Cjj 

(A-72) 

Ci ~ ^ ? wt cos 01 “ ^n,wt s t n 0L 

(A-73) 

Cm - C m 

(A-74) 

c n = c £,wt sin <* + c n,wt cos <* 

(A-75) 

C X a = - C D a cos 01 + sin a + °N = -C Aq , 

(A-76) 

C Xd = - C D d cos “ + C La sin a = -C A& 

(A-77) 
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C ?(i 

= % 






(A-78) 

Cy p 

“ Cy /3 






(A -79) 

c z a 

- _n T 
C D a 

sin oi 

- Ct cos a - 
• Lj q? 

C A = “ C N a 


(A -80) 

C Z & 

= -C T 

U CL 

sin a 

- Ct, . cos a = 

-Cn • 

iN 0? 


(A -81) 

%- 


cos a 

Cn /3,wt ' 

sin a 



(A-82) 

c ‘r 

c l- 

P,wt 

cos 01 ■ 

- C D(§;Wt sin a 



(A-83) 

%- 

C ^p ,wt 

cos^ a 

+ c n 

n r,wt 

sin^ 

a ” ( Cn P,wt 

+ Ci Y ^sin a cos a 

(A -84) 

c h = 

^r,wt 

0 

cos* a 

“ Cn p,wt 

sin^ 

a + (%,wt 

- C n ,Wn a cos a 
W r,wt/ 

(A-85) 

n 

O) 

II 

^6,wt 

cos a - 

■ C ns , sin a 
u 5,wt 



(A-86) 

Cm a 

= Cm fl 






(A-87) 

c m . 

CK 

= C m . 
a 






(A-88) 

Cm q : 

= Cm q 






(A-89) 

Cm 5 : 

= Cm 6 






(A-90) 


Cn /3 = %,wt Sin a + %,wt cos “ (A-91) 

S = %,wt Sin “ + S,wt C0S “ (A- 92 ) 

Cn P " %,wt cos2 a ‘ % >wt sin2 “ + (%,wt - C "r,wt) Sin “ cos “ (A- 93 ) 

c n r = C n r>wt cos 2 a + C ip wt sin 2 a + (Q r>wt + sin a cos a (A-94) 

° n 6 = %,wt Sin “ + Cn 6,wt C0S “ (A - 95 > 
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TRANSFORMATION EQUATIONS FOR MOMENTS AND PRODUCTS OF INERTIA 


Body to Flight Stability Axes 

i x s “ l x cos2 + l z sin2 «o " 2I xz sin "0 COS “0 

I Y,s = l Y 

l Z,s = l x sin2 <*0 + X Z cos2 “0 + 2I XZ sin a 0 cos a o 

*XZ,s = ( X X - I z) sin a 0 cos a 0 +I Xz( cos2 a 0 ‘ sin2 a o) 
I XY,s = *XY cos a 0 +I YZ sm a o 
*YZ,s = ^Z cos a 0 “ I XY sin a 0 

Body to Principal Axes 

X X P “ ix cos2 6 + ^Z s ^ n2 e ~ ^ixz s ^ n 6 cos 6 
I y? p = Iy 

I z p = I z cos 2 e + l X sin2 6 + 2I XZ sin 6 cos 6 


Flight Stability to Body Axes 
IX = *x,s cos2 a 0 + X Z,s sXn2 “0 + 2X XZ,s sin « 0 COS "0 


I Y = I 


Y,s 


I Z = ^X,s sin2 “0 + I Z,s cos2 “0 " 2I XZ,s sin “o C0S “o 

l XZ = (*Z,s ' I X,s) sin “o cos “o + I XZ,s( cos2 “0 " si " 2 «o) 
I Y z “ *yz,s C0S a 0 + X XY,s sin <*0 
X XY = I XY,s COS *0 - X YZ,s sin a 0 

Flight Stability to Principal Axes 

I Xj P = I X,s cos2 ” + I Z,s sin2 ’J + 2I XZ,s sin 71 cos 71 
l Y,P = 

I Z,P = X X,s sin2 71 + l Z,s cos2 V - 2I XZ,s sin * cos 11 


(A-96) 
(A -97) 
(A-98) 

(A-99) 

(A-100) 

(A-101) 


(A-102) 

(A-103) 

(A-104) 


(A-105) 
(A -106) 
(A -107) 

(A-108) 
(A -109) 
(A-110) 

(A -111) 
(A-112) 
(A-113) 
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Principal to Body Axes 


he = hc.p cos2 e + l z,p sin2 6 

(A-114) 

hr = hr,p 

(A-115) 

Ij = Ijj p sin^ e +Ig p cos 2 e 

(A-116) 

Ixz = (iz,p - hc.p) sin 6 cos 6 

(A-117) 


Principal to Flight Stability Axes 

!x,s = J X,P cos2 77 + I Z,P sin2 77 

l Y,s = ^P 

I Z,s = l X,V sin2 77 + I Z,P cos2 77 
X XZ,s = ( X X,P “ *Z,p) sin 77 cos 77 

SPECIAL FORMS OF EQUATIONS OF MOTION 

Coupled Linear Equations of Motion 

These are the linearized equations that describe the small perturbation motions of 
a vehicle about a steady -state flight condition. In steady -state flight, the components of 
thrust, aerodynamic force, and other forces in a given direction are balanced by the com- 
ponent of vehicle weight and, hence, the initial components do not appear in the equations. 
Other specific assumptions made in deriving these equations are 

(X) Total -velocity components, angular rates, and angles are equal to a steady - 
state value plus a small perturbation value; for example, the velocity com- 
ponent along the X-axis is u = U 0 + Au, where u is the total -velocity 
component, U 0 is the initial steady -state velocity, and Au is the per- 
turbation velocity 

(2) Products and squares of perturbations can be neglected 

(3) The vehicle has a plane of symmetry (Ixy ” lyz ” 9) 

(4) Pitch and roll perturbation rates are given in terms of the initial pitch and roll 

angles B Q and <p Q and not in terms of the total pitch and roll angles 
0 o + A6 and <p 0 + A <p; that is, 


(A-118) 
(A -119) 

(A-120) 

(A-121) 
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6 = Aq cos (f) Q - Ar sin 0 O 
and 

sin 6 0 sin 0 Q 

d> = Ap + Aq sin cb n + Ar cos cp 0 

F H ° cos 6 0 ° cos 6» 0 

(5) Aerodynamic cross -coupling coefficients are zero 

(6) Mass -damping and jet-reaction forces and moments can be neglected 

(7) The vehicle is flying at constant thrust 

These equations are written about the body axes but can be converted to the flight 
stability axes by letting W Q = 0 and Q 0 = y 0 . Also, to reduce the number of variables 
to six, the incremental velocity components Au, Av, and Aw can be expressed in 
terms of resultant velocity and incremental angle of attack by using the following equa- 
tions (the alternative would be to expand the aerodynamic forces on the right-hand side 
in terms of u,v,w derivatives): 

For body axis, 

Au = AV cos a Q - Ao? sin 

Av = Voo Aq? 

Aw = AV sin o?q + V^ A a cos o?q 
F or flight -stability axis, 

Au = AV 
Av = Voo A a? 

Aw = Voo A o? 



The six equations are 


u + Q 0 Aw + W 0 Aq - R 0 Av - V Q Ar + g ( J Aq dt cos Q q cos </) 0 - J Ar dt cos 0 O sin 




Cx V ^ + % fvi + 1 Cx 6 A5 + • • •) (A ' 122) 


Av + R 0 Au + Uo Ar - P 0 Aw ~ w o A P - g( Ap dt cos 6 0 cos <p 0 + J Ar dt sin 6> 0 j 

= a/3 + % ife + Cy p + ° Yr 

+ ^ Cy§ A6 + . . 


A6 + . . « 
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Aw + P G Av + V Q Ap - Q q Au - U Q Aq + g Ap dt cos 0 Q sin 0o + j Aq dt sin 0 O ^) 


pV °° S/ C z Ao' + C z .^l + C z ^X+C 7 

1 Zrv Z a 2Voo Z V Vco Z q 2Voo 


2m 


AS + . . . 


(A -124) 


Ap I x + (Qo Ar + R q Aq)(l z - I Y ) - (P Q Aq + Q 0 Ap + Ar)l xz 

-¥ v ^^ A ^ + %§z + %§z +c hwi + l c i6 AG + - • •) 

Aq ly + ^P G Ar + R 0 Ap^I x - I z ) + 2(P 0 Ap - R 0 Ar)l xz 
= | P VIS£ (c ma Aa + C mi + C ray + C m q + ]> C m g A6 + 


(A-125) 


(A -126) 


Ar I z + (P Q Aq + Q q Ap^Iy - I x ) + (Q q Ar + R 0 Aq - Ap)l xz 


-ipvfWc^ A£ + C n - Ml 


Ap l 


$ 2V X + Cn P 2V qo t ^ n r 2V m 


+ C 


Ar 1 


i 


+ / c n 5 AS + . . . 


(A-127) 


Uncoupled Equations of Motion 

The linearized equations of motion can be uncoupled (lateral motions made inde- 
pendent of longitudinal motions) by assuming that the vehicle is in straight and level flight 
and that there are no components of initial velocity except U Q , W 0 , and Q 0 in the 
initial steady-state condition (i.e., V 0 = P 0 = R 0 = Tp Q = <fi Q ~ 0). Under these assump- 
tions, longitudinal equations contain only the variables Au, A 0, and Aw (q = 0); and 
lateral equations contain only the variables Av, Ar, and Ap. 


Longitudinal equations . - 

2 / 

Au + Qo Aw + W 0 Aq + g% cos e g;0 = ( C Xq Aa + C Xa + C Xl AV 


2V V< 




AS 


(A -128) 


78 



APPENDIX A - Concluded 


pvl s 


Aw - Q 0 Au - U 0 Aq + g <9 g sin 6 g>0 = Aa + C z ^ 


AafL 


AV 


+ C7 - 

a 2Voo Z V Voo 


+ C Zq^f + 2 Cz 6 A5 , 


(A-129) 


Aq I Y = £ ^(Cm Q , Aot + C raA 


Lateral equations. - 


A al , r, AV , r, Aq £ 

• "oTr — + 77 — + 0 ,r“ 

o? 2V 00 V Voo Q 2VOO 


+ / Cm* A6 


Av + U 0 Ar - W 0 Ap - g^0 g cos 0 g>o + xf/ g sin 0 gjO ^ 

” ( Cy /3 ^ + Cy |3 + % fvf + C Yr WZ + I ' ° Y S A6 ) 


pv^s, 

2m 1 


(A -130) 


(A-131) 


Ap lx + Q 0 Ar ( X Z - *y) * (Qo Ap + Ar)l X 2 

pV^Si ^ A/3 £ , r< Ap i , ^ 

= — 2— t Ci )3 A(3 + %2VZ + % 2V^ + Ci 


Ar * + Y<* fi A8 


r 2V 


(A-132) 


r I z + Q 0 A p(!y " *x) + (Qo Ar “ A p)^XZ 

■ 2 < 

oo* 

2 


pV^Slfr, A „ A/3 J t p Api , r, Art Vf a) 

= — i~ ( Cn /3 A(3 + Cn /3 2V^ ° n P + Cn r 2V^ + Z Cn 5 A5 y 


(A-133) 


Wind-Axes Equations for a Point Mass 

The wind -axes equations used in ballistic trajectory studies in which the vehicle is 
considered to be a point mass are, along the flight path, 


m — = C D - W sin y + T 

dt 2 li 


(A -134) 


and, normal to the flight path, 


mV^= C L - W cos y 

dt 2 


(A-135) 


These equations assume a constant thrust acting parallel to the flight path and 
neglect jet damping and reaction control forces. 
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DERIVATION OF EQUATIONS PRESENTED IN SECTION II 
FOR TRANSFER OF COEFFICIENTS TO NEW 
REFERENCE CENTER 


The aerodynamic forces and moments acting on a flight vehicle can be considered 
to be functions of six independent variables: a ,/3,V,p,q,r in one system, or u,v,w,p,q,r 

in the alternate system. To evaluate force and moment derivatives at a new reference 
center, such as at a new c.g., the derivatives at the new center, which represent changes 
in forces or moments with changes in the independent variables, both evaluated at the new 
center, are written, for example, 


ax’ _ 

_ ax’ 

da 


ax’ a/3 . 

ax’ 

av 

ax’ 

ap 

+ ax 

aq 

+ 9 X 

T 

3r 

da' 

da 

da 1 


3/3 da' 

av 

da' 

ap 

da' 

aq 

da' 

ar 

da 

ax „ 

_ 3X. 

da 

+ 


ax’ 


ax’ 


ax’ 


ax' 

dr 

a/3 1 

da 

a/3 1 


a/3 a/3 T 

av 

dp' 

ap 

3/3 

aq 

dp' 

ar 

dp' 

ax’ _ 

_ 9X 

da 

+ 

&LM. + . 

ax' 

^ + 

ax' 

i£. + 

ax’ 


ax' 

dr 

ar 1 

da 

dr' 


a/3 ar T 

av 

ar 

ap 

3r 

aq 

dr 

3r 

dr' 


where the forces X ,Y ,Z and moments L ,M ,N at the new reference center are 
related to those at the original reference center by 




X =x 

L = L + zY - yZ 

Y = Y 

M’ = M + xZ - zX > 

z’ = z 

n’ = N + yX - xY j 


The variables a',/3,V,p,q,r evaluated at the original reference center can now be 
expressed as dependent functions of the independent variables a ,p' ,V T ,p’ ,q’ ,r’ of the 
new system through the equations 
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0 ! = tan-1 2 = tan- 1 V|, sin ^ cos p' -py + qx 
u V cos a cos p - qz + ry 


p = sin’ 1 'L = sin’ 1 v 

V V 


> (B*-3) 


V T = Ju 2 + v 2 + w 2 

= ^j(v' cos at’ cos p' - qz + ry)^ + (v' sin p' rx + py^ + (v ? sin a' cos /3 f - py + qx)^ 



(B-4) 


Since from equations (B-4) and derivatives such as 

9p f dp’ 9q’ 8q’ 9r T 8r* 

= 0 (they are derivatives of one independent variable with respect to another), 
the derivatives in equations (B-l) can be evaluated from equations (B-3) and (B-4) as 


dot 

da' 

Jg. 

9/3’ 

9V 

9V f 

9a; 

9p’ 

9p_ 

9P’ 

dot 

9q* 

J*L 

9q’ 


= 1 


= 1 


= 1 


-y cos a 
V cos /3 


= 1 


x cos a + z sm a 
V cos /3 


= 1 


da -y sin 
9r* " 


V cos /3 


9r 

9r T 


= 1 


9V 


9p 


- = z sin /3 - y sin a cos j3 


9/3 = _9V_ _ _9p_ _ _9q_ _ _9r_ = Q 
da' da' da' da' da' 

da _ dy _ dp _ dq _ dr _ Q 
9/3* ~ dp' ~ dp' ~ dp' "* 9/3* ~ 

90! _ dp _ dp _ 9q _ 9r _ n 
9V’ *" 9V’ ” 9V’ ~ 9V’ ~ 9V’ ” 

dp _ z cos p + y sin a sin p 

dp' y 

= J* = o 

9p f 9p T 


M. - (*. °f - y.an.,«)sln,J i®5L = (x sin a - z cos a) cos 
9q V ,a q’ 

9p _ 9r - 0 

9q’ 9q’ 

M. = -X COS P - y cos a sin p SV_ = - „ ^ _ - sJn 0 

Sr* 1 y 17 r- 


> (B-5) 


9r 

_3e_ = JSL = 0 

9r’ 9r’ 


9r’ 
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Substituting equations (B-5) and (B-2) into equations (B-l) and nondimensionalizing yields 
the transformation equations given in section n. 

The transformations for the u,v,w derivatives (eqs. (11-41) to (11-77) of Sec. n) 
are derived in the same way except that in this system the derivatives are of the form 


ax' = ax' au ax' av ax' aw ax^_ _ap_ ax_ _a^_ _ar_ 

au’ au au’ + av au’ + aw au’ + ap au’ aq au T ar au' 


> (B-6) 

ax' _ ax' au ax' ay ax' aw ax' ap t ax aq ( ax ar 

ar 1 - au ar’ + av ar’ + aw ar T + ap ar’ aq ar’ ar ar’ 


The derivatives in equations (B-6) are evaluated from equations expressing velocity com- 
ponents at the original reference center in terms of those at the new reference center as 

u = u’ - qz + ry 

v = v’ - rx + pz ) (B-7) 

w = w T - py + qx 
and are given as 


j 


au 

av 

aw 

_ 9p 

= J£L = . 

au T 

” av’ ' 

" aw’ 

ap’ 

aq’ 1 

au 


av 


aw - 

ap’ 

= 0 

ap’ 

= z 

8p' - " y 

au 

= -z 

av 

= 0 

^ = 5 

aq’ 


aq’ 


aq* 

au 

= V 

av 

= -X 


ar’ 

J 

ar’ 


ar’ 


(B-8) 


All other derivatives are zero. 

Substituting these derivatives from equations (B-8) along with equations (B-2) into 
equations of the form (B-6) and nondimensionalizing then yield the transformations for 
u,v,w derivatives given in section II. 

In nondimensionalizing, forces are divided by q^S; moments, by q^Sj^; p, q, r, 
a, and j3 are multiplied by 1/2 V and u,v,w by 1/V. The velocity derivatives, such 
as C Xv , Cy v , and C mv , are written in terms of static forces and moments by making 
use of the following relationship: 
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X = ipV 2 SC x 


therefore, 


H = pvsc x + ipV 2 S 


8C X 

av 


(B-9) 


The derivative dC^^V represents the change in Cx brought about by changes in 
velocity V by itself as an independent variable and not by changes in velocity due to 
angular -velocity components p,q,r for the rotating body. They are the result of such 
things as aeroelastic and Mach number effects and are usually neglected, so that equa- 
tion (B-9) becomes 

|f=PVSC x (B-10) 

The relationship in equation (B-10) is also used in nondimensionalizing the derivatives. 

The transformations given in section II for the static force and moment coefficients 
and the static -stability derivatives (derivatives with respect to a or p) are simplified 
forms that apply only if there is no vehicle rotation and p, q, and r are zero. Although 
it probably is not practical to use such forms, more general forms can be derived that 
account for significant rotation; for example, it can be shown that the complete form of 
da /da 1 from equations (B-3) is 


= 1 + (qz - ry) cos a + (py - qx) sin a 
da 1 

and, by similarly evaluating dp/ da 1 and 3V/8 q! t , it can be shown that the more com- 
plete transformation for Cxq, is 


c x a = c x a 

- c Xj3 
+ 2C X 


l + (qz - ry) cos a + (py - qx) sin a 


V cos p 

(qz - ry) sin a - (py - qx) cos a 
V 

(qz - ry) sin a - (py - qx) cos a 


V 


sin p 
cos p 


(B-ll) 


which reduces to = Cx (the form given in Sec. II) if p = q = r = 0. 
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METHODS OF MEASURING CENTER -OF -GRAVITY LOCATIONS AND 
MOMENTS OF INERTIA OF MODELS AND FLIGHT VEHICLES 

The methods presented here are a summary of methods used at the Langley 
Research Center to measure moments of inertia of rocket -propelled models and the 
methods presented in reference 3 for determining c.g. locations and moments of inertia 
of full-scale airplanes. In these methods all oscillations are assumed to be small. 

Some other limitations and precautions that must be taken in using these methods, and 
which apply to all methods, are discussed in reference 3. 

CENTER -OF -GRAVITY LOCATION 
Longitudinal c.g. Location 

The longitudinal c.g. location can be determined by mounting the aircraft or rocket 
on weighing scales. For a typical aircraft installation, two scales measuring the forces 
R 1 and R2 are located at wing jack points and a third scale measuring the force R3 
is located at some distance l forward or aft of the jack points. The longitudinal dis- 
tance from the jack points to the c.g. is then given by 

R3I 

X ~ Rj + R2 + R 3 ( C-1 ) 


Vertical c.g. Location 

In order to determine the vertical c.g. location the vehicle is mounted on a knife 
edge as shown in figure 10. The weighing cradle has a weight W c and a centroid z c , 
and the vehicle is supported in various roll and pitch attitudes by a vertical reaction 
force Rj acting at a distance y^. The distance from the knife edge to the airplane c.g, 
is obtained from measurements of R^, y^, and 0 by using the equation 

R lYl - W c z c sin 0 


z = 


W sin 0 


(C-2) 


MOMENTS OF INERTIA 


Compound -Pendulum Method (Fig. 11) 

This method is used mainly to determine pitch and roll moments of inertia on small 
models that can easily be mounted with a single attachment point. The moment of inertia 
of the model about an axis through the model c.g. is given by 
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T W M+S z M+S T) 2 W S z S p 2 W M Z M , r _ s \ 

• S ' * 

The subscript M+S refers to the model plus the supporting hardware, S to the sup- 
porting hardware alone, and M to the model alone. The distance z ]y[.^g’ from the knife 
edge to the c.g. of the model plus supporting hardware, is obtained from 


Wg z s + w m^M 


J M+S 


W 


M+S 


(C-4) 


Spring Method 

This method is usually used to measure pitch and roll moments of inertia when the 
model is so large that more than one support is needed. The model is mounted on a knife 
edge and springs are attached at equal distances on both sides of the knife edge, as shown 
in figure 12. The moment of inertia can be determined from the equation 

I = Mi p2 + E££ P 2 . m5 2 (C-5) 

477-2 47J-2 


Spring Method for Full-Scale Vehicles 


This method, which is described in detail in reference 3, uses two sets of springs 
with each set having a different spring constant, so that both the moments of inertia and 
the vertical c.g. location can be determined. The springs are arranged as shown in fig- 
ure 13. The moment of inertia about the point axis through the knife edge is given by 


I = 


Ft.i - K t, 2 ) p j 



-2 

- mz 


- Ic 


- m c z 


2 

c 


(C-6) 


The vertical distance from the knife edge to the c.g. of the model is then determined from 



(C-7) 


where W c and z c are the weight and c.g. vertical displacement (measured from the 
knife edge) of the weighing cradle that is used. The constants Kj. ^ and are 

determined from the constants for the springs shown in figure 13 by the equations 
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*^,1 = 2K s,l a s 2 
^,2 = 2K s,2 a s 2 


Torsion-Pendulum Method 


Used for measuring yaw moments of inertia, this method involves mounting the 
model on a torsion shaft as shown in figure 14(a). The moment of inertia about the lon- 
gitudinal axis of the shaft is given by 


I = 


kP 


M+S 


47T 2 


kp£ 

47T 2 


(C-8) 


Multifilar -Pendulum Methods 

Essentially the same as the torsion spring method except that the vertical wires or 
rods provide the restraining spring moment, an example of a bifilar pendulum is shown in 
figure 14(b). The moments of inertia are determined from 

j _ W M+S P M+S j2 W S P S 2 -2 (c 

16*2 £’ 1 ^ 2 ~ & (C ' 9) 

where a is the diameter of the circle around which the wires or rods are attached and 
is the length of one of the wires or rods. A system with two wires is called a bifilar 
pendulum; a system with three is called a trifilar pendulum. In using these methods, the 
wires or rods must be centered about the system mass center. 
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DETERMINATION OF LONGITUDINAL AND LATERAL STABILITY 
DERIVATIVES BY USING SIMPLIFIED LINEAR ANALYSIS 

4 


Perhaps the simplest and most direct method of analyzing transient flight data is to 
determine static- and dynamic -stability derivatives from the frequency and the time to 
damp to half -amplitude of transient motions, such as the damped angle -of -attack oscilla- 
tion shown in figure 15. This method has been used effectively to determine the control 
response characteristics of rocket -propelled models in reference 9. Even when a more 
sophisticated analysis is planned, this method can provide, before the detailed analysis 
is begun, a fairly accurate quick-look assessment of the vehicle characteristics from the 
basic flight records. 

The oscillation in figure 15 represents the transient response to a step -control 
input. The method assumes that 


(1) The forward velocity and the Mach number of the vehicle are constant 


(2) The longitudinal aerodynamic forces vary linearly with angle of attack a, pitch 
angle 0, elevator deflection S e , dce/dt, and d0/dt; lateral aerodynamic forces are 
assumed to vary linearly with sideslip angle /3, yaw angle \ //, rudder deflection 6 r , 
d/3 /dt , and d^/dt 


(3) The vehicle is in level flight before the control deflection is applied 

(4) Longitudinal and lateral motions are independent of each other (all aerodynamic 
cross derivatives such as Q q and C mp are neglected) 

Under these assumptions, the longitudinal equations of motion can be simplified to 


mV oo/d0 dee 


<1 0C S V dt 


^ - c T n + Ct ol + C L 

dt dt L >° u ot u 


6 


(D-l) 


X Y d 2 g 

1oo S5 dt 2 


= C 


m,o 


+ C m ol + C 


a c 


OL 


m 


a 2V C 


+ Cm 


qc 


<1 2V r 


+ Cmg 5e 


(D-2) 


where Cl 0 and ^m,o are va -l ues of lift and pitching -moment coefficient in 
trimmed level flight. The solution to equations (D-l) and (D-2) is of the form 

ol = Ce a ^ cos (cut - £) + (D-3) 

In equation (D-3), C is a constant that is determined from initial conditions, ce^ is the 
trim angle of attack, 


87 



APPENDIX D — Continued 


a =-2 


c L 


m 


“ - ( Cm q + Cm a)Wc J 


(D-4) 


'111 


and 


w "\ 

| X1J 

m T = 


q» s 

V = - 



a 


C m C I 

5 q “ - (a-) 2 


2V,*, I’m’ 


' (D-5) 


(D-6) 


q.ooSc 


The damping constant a' and the frequency c o are determined from measured quanti- 
ties; a’ can be calculated from the rate of decay of the oscillation (see fig. 15) as 


a* = 


A (*2 

l0g e lo S e I -0.693 


(D-7) 


*2 - *1 * 1/2 * 1/2 
is the time for the oscillation to damp to one -half its initial amplitude. The 


where t-jyg 

frequency c o can be calculated from the period of the oscillation P as 
2i t 


co = 


Equation (D-5) can be rewritten as 


C - I’ffa') 2 , V 2 1 £ CmqCLa 

C m „- 1 L' a > + a, J 2 V M m' 


(D-8) 


(D-9) 


and equation (D-4), as 


4I T V, 


c L 


Cmq + C m . r — la* + 


oi 


OL 


2m 1 


(D-10) 


The last term of equation (D-9) is usually very small compared to the first term (usually 
less than 1 percent) and may be omitted in most cases. (See ref. 9.) Omitting this term 
and substituting equations (D-7) and (D-8) give the equations for the longitudinal static - 
and dynamic -stability derivatives in terms of the period and the time to damp to half- 
amplitude of the oscillation as 

,2 


'm 


a 


a 2 t 
- 47 r Iy 

1 

q Sc 1 

°° I 

p2 


1 /0.693^ 


p2 ^Tr^yl/2 


(D-ll) 
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-2Iy fl.386V cq 

(c mq + c m . ) = 



(D-12) 


Similar equations for the lateral derivatives can be developed as 

,2 


Cr» ^ — 


4j t 2 1 z 

1“ 

1 

q Sb 

oo 

f 2 

- C n . N 
n i 3y 



1 /0.693 


P^ 47T^^i/2 


-21 


z / i.386V„ g 
: l qJ3t. w Y 


1/2 


(D-13) 


(D-14) 


For the lateral-derivative equations a small angle of attack is assumed. The values 
of C n p and ^Cn r - C n ^ given by equations (D-13) and (D-14) are those about the stabil- 
ity axes. 

It should be noted that the mean values of angle of attack or sideslip do not neces- 
sarily correspond to the trim value, so that the mean value that is used to determine 
ti/g (eq. (D-7)) must be obtained by selecting positive and negative amplitudes from the 
angle -of -attack or sideslip envelopes and determining the mean of these positive and 
negative amplitudes. 
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USE OF DIRECTION COSINES AND QUATERNIONS 
IN MOTION CALCULATIONS 

In the equations of motion and throughout this report, vehicle orientation with 
respect to given axes is specified in terms of the Euler angles Two other 

methods can be used, however, to specify alinement. One is to give the alinement in 
terms of direction cosines; the other is to specify alinement in terms of ’’Euler param- 
eters,” which are components of a four -parameter quantity called a quaternion. Whereas 
there is no particular advantage in one of these methods over the other, both have certain 
advantages over the use of the Euler angles. Both eliminate the singularities that occur 
when vehicle attitudes approach ±90° (a condition known as gimbal lock) and, whereas the 
equations relating Euler angles to angular rates are nonlinear (eqs. (V-7) to (V-9) of 
Sec. V), those relating direction cosines and Euler parameters to angular rates are linear, 
so that the computational procedure is simplified, particularly in analog computations. 
Either method is preferred to the Euler angle method, therefore, in certain applications. 
The basic equations needed in applying these two methods are given here. 


DIRECTION -COSINE METHOD 

In the direction-cosine method, the axes transformation for any vector has the form 

(E — 1) 


"x ,_ 


X" 

Y’ 

= D 

Y 

_z’_ 


Z_ 


where D is the matrix made up of direction cosines and is defined as 


D = 


where the elements of D are given in terms of Euler angles as 


*xx 

l XY 

l XZ 

£ yx 

%Y 

l YZ 

^zx 

l ZY 

Q-zz 


(E-2) 
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^XX " cos ^ C0S ® 

^XY = s ^ n ^ COS ^ 
j^XZ “ _s i n 0 

£ yx = -sin xp cos ^ + cos \p sin 6 sin 0 
j^ YY - cos i// cos <jf) + sin \p sin Q sin p ^ 
£ Y z = cos 0 s ^ n 0 

^ZX “ s ^ n ^ s ^ n ^ + cos ^ s ^ n ^ C0S ^ 
^ZY ” ~ cos ^ sin 0 + s ^ n ^ s ^ n 0 cos 0 
#ZZ — cos 0 cos (p 


(E-3) 


It can be shown that D can be determined from the angular velocities by using equa- 
tions (E-4) and (E-5) as follows: 


0 

1 

CO t 

z t 

-®r 

^z' 

0 

“x’ 

V 

_c v 

O 

I 


which is a linear (matrix) differential equation. The angular velocities 


(E-4) 


co 


X 


.00 



are those of the primed axes system (see eq. (E-l)) with respect to the unprimed axes 
system. In motion calculations where the primed axes are considered as the vehicle 
body or other reference axes and the unprimed axes are taken as the gravity-axis sys- 
tem (see fig. 5), oi „a> ,,w , can be determined from the body-axes angular velocities 

P,q,r by 


00 , 


P~ 


-L 

X 





V 

= 

q 

- D 

-(p e + X) cos L 

00 t 

L z_ 


r 


-(pe + x) sin L 


(E-5) 


or, for a flat nonrotating earth, 




v = p 



(E-6) 
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QUATERNION METHOD 


Basic Quaternion Relationships 

In the quaternion method an axes transformation has the form (ref. 10) 


x T 


— i 

X 

Y’ 

= G 

Y 

z'_ 


A 


(E-7) 

where G is a matrix composed of the Euler parameters e o? e i> e 2? e 3 and is defined as 
( e 0 2 + e l 2 - e 2 2 - e 3 2 ) 2 ( e l e 2 + e 0 e 3) 2 ( e l e 3 - e 0 e 2) 

2 ( e l e 2 ' e 0 e 3) ( e 0 2 - e l 2 + e 2 2 ' e 3 2 ) 2 ( e 0 e l + e 2 e 3) I ( E ' 8 ) 

2 ( e l e 3 + e 0 e 2) 2 ( e 2 e 3 ‘ e 0 e l) ( e 0 2 ' e l 2 " e 2 2 + e 3 2 ' 

The Euler parameters are elements of a four -parameter quantity called a quaternion, 
defined as 


G = 


q = e 0 + e^ + e 2 j + e 3 k 

where e 0 ’ e l’ e 2’ e 3 are real numbers and the vectors i,j,k satisfy the following 
conditions: 

i 2 = j 2 = k 2 = -1 
ij = -ji = k 


(E-9) 



(E-10) 


The quantity e q is the real or scalar part of the quaternion; the terms e^i + e 2 j + e^k 
make up the imaginary part. The length or norm of the quaternion is defined as 

| q | “ = \j e o 2 + e i 2 + e 2 2 + e 3 2 (E -11) 

The quantity q* is the conjugate of the quaternion and is defined as 

q* “ e 0 - e i^ " e 2 3 " e 3 k (E-12) 

The transformation in equation (E-7) can also be written in terms of the quaternion and 
its conjugate as 


(E-13) 


x r 


~X~ 

Y T 

= q* 

Y 

_Z T _ 


Z 
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In flight -motion calculations, in which the Euler parameters define the alinement 
of the vehicle body or other reference axes with respect to the gravity axes (see figs. 4 
and 5), the differential equation from which the Euler parameters can be determined is 


“«0 


" e l " e 2 " e 3 


r 

p~ 


-L 

a 

e l 

e 2 

ii 

COll- 1 

e 0 " e 3 e 2 

e 3 e 0 ' e l 

< 


q 

r 

- G 

-{p e + X.) cos L 
+ Xjsin L 

> 

„ e 3_ 


“ e 2 e l G 0_ 


V 



\ ^ / __ 



4 

(E-14) 


where G is defined by equation (E-8). 


Relationship Between Euler Parameters and Euler Angles 


Initial values of e o> e i> e 2’ e 3 for use in P erformin & the calculation indicated by 
equation (E-14) can be determined from the initial values of the Euler angles. The Euler 
parameters are related to Euler angles by the following equations: 




t 


eg = cos ~ cos ^ cos 4 + sin 4 sin ^ sin ^ 

= cos cos | sin £ - sin sin | cos |- 
e 2 = cos ^ sin | cos |- + sin & cos | sin |- 


e 3 = sin cos | cos ^ - 


COS 




a 

sin ~ sin 


4 


(E-15) 


sin 6 - 
tan 0 = 


tan 0 = 


~ 2 ( e l e 3 ' e 0 e 2) 

2 ( e 0 e l + e 2 e 3) 

p2 +e 2- e 2-e2l 

e 0 1 2 3 > 

2 ( e l e 2 + e 0 e 3) 

p 2 „ e ^ - e ^ + e« 2 
e 0 e l e 2 +e 3 J 


(E-16) 
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APPENDIX F 


SCALING PARAMETERS 


This appendix summarizes parameters for the three types of scaling used in aero- 
dynamic testing: dynamic, aerodynamic, and aeroelastic. These parameters establish 
the requirements for a given type of scaling in that each significant parameter must have 
the same value for both model and prototype if the simulation is to be valid. In most 
cases, it is impossible to satisfy all the scaling requirements at the same time and com- 
promise techniques have to be worked out. A knowledge of the scaling parameters is 
essential to understanding and applying these techniques. 

DYNAMIC -SCALING PARAMETERS 

Dynamic -scaling parameters established the conditions under which the motions 
(accelerations, velocities, and angles) and the forces for the model simulate those for 
the prototype. These parameters are obtained by requiring that, in the equation describ- 
ing the motion of the flight vehicle, the ratio of any one term to another has the same 
value for model and prototype. The equation of motion has the general form 


Inertial force = Weight + Aerodynamic force 


(F-l) 

Two scaling parameters determined from equation (F-l) 

are 


_ Vehicle inertial force 

mv|/e 

V? 

/ TTl <)\ 

Vehicle weight 

w 

gl 

{V -Z) 

Mass ratio - Y e ^ c ^ e inertial force 
Aerodynamic force 

mvJ/4 

pV 2 l 2 

w 

pg£ 3 

(F-3) 


In addition, the scaling assumes geometric similarity between model and prototype 
as well as similarity in mass distribution. The scale factors given in table VI were 
obtained by satisfying all these requirements. 


AERODYNAMIC -SCALING PARAMETERS 

Aerodynamic -scaling parameters must have the same value for model and proto- 
type if the flow field around the model and, hence, the aerodynamic force and moment 
coefficients, is to be the same as around the prototype. Aerodynamic -scaling param- 
eters determined from the equation of motion of the fluid, in which viscous, pressure, 
and gravity forces are considered, are 
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APPENDIX F - Continued 


Reynolds No. 


Fluid inertial force 
Fluid viscous force 


Mach No. 


Fluid inertial force 
Fluid pressure force 


Froude No. = 


Fluid inertial force 
Fluid gravity force 


_ pO z pv c 0 <> 

fiVoof ' 4 

(pvyV 2 Xc 

[p^i* ) ~ a 

(F-4) 

(F-5) 

$ 

pv|f 2 _ v| 

pgJ^3 ^ 

(F-6) 


In addition, if surface-tension forces are important, a fourth parameter to be con- 
sidered is 


Weber No. 


Fluid surface-tension force _ ol _ o' 
Fluid inertial force pV<^£ 2 pV 2 £ 


(F-7) 


where o is the surface tension per unit length. 

The pressure force considered in Mach number (eq. (F-5)) is that due to the pres- 
sure differential across a shock wave in compressible flow. The pressure -force param- 
eter considered in incompressible flow is 


Euler No. 


Fluid pressure force _ P 
Fluid inertial force pVj 2 & 2 


(F-8) 


The Euler number is usually not important because usually the resultant body forces are 
measured; however, it becomes important when body forces are determined from mea- 
surements of pressure distribution. 


Two scaling parameters determined from the laws of thermodynamics are 


Prandtl No. 


1 Heat added by convection 

Reynolds No. Heat added by conduction 


v C P pV °° T A _ c p/ d (F-9) 

~ VJt kT/jt 2 k 

where v - p/p, k is the coefficient of heat conduction (thermal conductivity) of the 
fluid, and 

Grashof No. = (F-10) 

vjTo 

where AT is the temperature difference between two representative points in the fluid 
and T 0 is a representative temperature. 
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APPENDIX F - Continued 


Aeroelastic -Scaling Parameters 

The basic aeroelastic -scaling parameters to be satisfied in simulating deforma- 
tions caused by aerodynamic loads are 


For structural elongation: 
EA 
pvie 2 

For structural bending: 


El 

pV 2 f 2 

^ oo 


(F-ll) 


(F-12) 


in which it is assumed that the strains in the model structure are the same as for the 
prototype; that is, 


1“ = 1.0 (F-13) 

e P 

It is also assumed that neither model nor prototype materials are stressed beyond their 
elastic limits and that the ratios of shear -to -tensile stress and of shear -to -elastic 
modulus are the same for model and prototype. A more thorough discussion of elastic 
scaling is given in reference 11. 


COMMENTS ON SCALING PROCEDURES 


Dynamic Scaling 

In general, tests of dynamically scaled models are limited to subsonic speeds 
(below M = 0.6) because of the requirement that the force coefficients must be the same 
on model and prototype. For vehicles designed to fly at transonic or supersonic speeds, 
the model is usually tested at the correct full-scale Mach number. The motions, loads, 
and accelerations of the test models are then reduced to nondimensional coefficients by 
utilizing the known mass and inertial characteristics of the models, and these coeffi- 
cients are in turn used to calculate the motions, loads, and accelerations of the prototype. 
In any case, care must be taken to insure that Reynolds number effects do not introduce 
an unacceptable distortion. Reynolds number effects can be neglected in most tests if 
the model Reynolds number is above 10^; however, in general, the Reynolds number 
effects on skin-friction drag must be accounted for. 

Aerodynamic Scaling 

In wind-tunnel and free -flight model testing, it is impossible to satisfy all the 
aerodynamic -scaling parameters at the same time; however, the problem usually reduces 
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APPENDIX F — Concluded 


to one in which only two of the parameters, Reynolds number and Mach number, are 
important. In most problems, the effect of the fluid gravity on the motion is negligible, 
and so the Froude number and the Grashof number of the fluid can be neglected. It can 
be shown that the Prandtl number is the same for model and prototype if the ratio of 
specific heats cp^/c v is the same for both, and this condition is automatically satisfied 
if air is used as the fluid in the model tests. A still further simplification is possible 
because Reynolds number effects are limited to the narrow boundary -layer region at the 
body surface and the flow outside the boundary layer is practically without the influence 
of viscosity. Thus, Reynolds number and Mach number effects can be considered sepa- 
rately. Most wind-tunnel tests at transonic and supersonic speeds are made at the full- 
scale Mach number. The Reynolds number effects, called scale effects, are accounted 
for by making corrections to skin -friction drag coefficient (the coefficient primarily 
affected by Reynolds number) or by using boundary -layer trips, wires or a rough grit, to 
cause separation at the proper chordwise position and therefore to simulate the full-scale 
flow pattern. These techniques are discussed in textbooks on wind-tunnel testing (for 
example, ref. 12) and reports dealing with their application, such as reference 13. 
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TABLE I.- CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


Physical quantity 

Multiply 

U.S. Customary Unit 

by 

conversion factor 

to obtain 
SI Unit 



(a) 

(b) 

Angle 

deg 

0.01745329 

rad 

Length 

ft 

0.3048 

m 

Velocity 

ft/sec 

0.3048 

m/sec 


knot (international) 

0.5144444 

m/sec 


knot (U.S. statute) 

0.44704 

m/sec 


mph 

0.44704 

m/sec 

Acceleration 

ft/sec 2 

0.3048 

m/sec 2 

Pressure 

lbf/ft 2 

47.88026 

N/m 2 


atmosphere 

1.01325 X 105 

N/m 2 


ft of H 2 0 (39.2° F) 

2.98898 X 10 3 

N/m 2 


in. of Hg (600 F ) 

3.37685 X 10 3 

N/m 2 

Moment of inertia 

slug -ft 2 

1.355818 

kg-m 2 

Mass 

slugs 

14.59390 

kg 

Force 

lbf 

4.448222 

N 

Moment 

ft-lbf 

1.355818 

N-m 

Area 

ft 2 

0.09290304 

m 2 

Mass flow 

slugs/sec 

14.59390 

kg/sec 

Density 

slugs/ft 3 

515.3788 

kg/m 3 

Spring constant 

lbf/in. 

175.1268 

N/m 


lbf /ft 

14.59390 

N/m 

Gravitational constant 

ft 3 /sec 2 

0.02831685 

m 3 /sec 2 

Specific heat 

Btu(thermochemical) 
lbm -of 

4.184 X 10 3 

joule/kg-°C 

Kinematic viscosity 

ft 2 /sec 

0.09290304 

m 2 /sec 

Coefficient of 

Btu(thermochemical) -in. 

518.87315 

j/m-sec-K 

heat conduction 

ft 2 -sec-°F 


Temperature 

o F 

(c) 

K 


°C 

(c) 

K 


°r 

(c) 

K 

Coefficient of 

slugs/ft-sec 

47.880258 

N-sec/m 2 

viscosity 





a Based on values in ref. 14. 

b Prefixes to indicate multiples of SI Units are as follows: 


Prefix (symbol) 

Multiple 

kilo (k) 

10 3 

hecto (h) 

10 2 

deka (da) 

10 

deci (d) 

o 

1 

h-* 

centi (c) 

lO- 2 

mill! (m) 

10~ 3 

micro (jli) 

10 " 6 


c Temperatures related by following formulas: 
K = |(OF + 459.67) 

K = OC + 273.15 
fC-|°R 
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TABLE III.- ANGLE DESIGNATIONS FOR DIRECT a TRANSFORMATIONS 


Axes conversion 

Angle designation 

From 

(Unprimed coefficients) 

To 

(Primed coefficients) 

A 

B 

Principal 

Body 

e b 

0 

Flight stability 

Body 

OiQ 

0 

Flight stability 

Principal 

^b 

0 

Wind-tunnel stability 

Body 

a 

0 

Wind-tunnel stability 

Principal 

a - e 

0 

Wind-tunnel stability 

Flight stability 

a - 

0 

Wind 

Body 

a 

0 

Wind 

Principal 

a - e 

0 

Wind 

Flight stability 

a - a 0 

0 

Wind 

Wind-tunnel stability 

0 

0 


a Direct transformations (eqs. (Ill -I) to (III-54)) represent rotation from ref- 
erence (unprimed) axes system through angle -B about Z-axis and then through 
angle +A about Y-axis to new (primed) axes system. Inverse transformations 
(table IV and eqs. (III-55) to (HI-107)) represent reverse of direct. 

b Angles e and r) between the principal axes and the body and flight stabil- 
ity axes, respectively, can be determined from moments of inertia by 


tan 2e 


tan 2?7 


2 ixz 
Iz - lx 

21 XZ,s 
^XjS ~ *Z,s 
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TABLE IV.- ANGLE DESIGNATIONS FOR INVERSE a TRANSFORMATIONS 


Axes conversion 

Angle designation 

— — i — 

From 

(Primed coefficients) 

To 

(Unprimed coefficients) 

A 

B 

Body 

Principal 

e 

0 

Body 

Flight stability 

a 0 

0 

Body 

Wind-tunnel stability 

a 

0 

Body 

Wind 

a 

0 

Principal 

Flight stability 

T] D 

V) 

0 

Principal 

Wind-tunnel stability 

o' - e° 

0 

Principal 

Wind 

a - e 

0 

Flight stability 

Wind-tunnel stability 

a - a q 

0 

Flight stability 

Wind 

a - ot Q 

0 

Wind-tunnel stability 

Wind 

0 

0 


a Inverse transformations (eqs. (HI-55) to (HI-107)) are reverse of direct 
(table III and eqs. (IH-1) to (HI- 54)) and represent rotation of primed axes system, 
through angle -A about Y-axis and then through angle +B about Z-axis, until 


axes coincide with original unprimed system. 

k Equations for determining angles e and rj given in footnote in table III. 
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TABLE V.- DESIGNATIONS OF FORCE AND MOMENT COEFFICIENTS 
FOR DIFFERENT AXES SYSTEMS 


Component 

Coefficients for axes system — 

Body or 
principal 
(a) 

Flight 

stability 

Wind-tunnel 

stability 

Wind 

X-axis force 

C x or -C A 

c X,s 

-C ? 

- C D 

Y-axis force 

C Y 

C Y,s 

c Y 

C C 

Z-axis force 

C z or -C N 

C Z,s 

' C L 

-c L 

X-axis moment (roll) 

C l 

C l ,s 

^ l ,wt 

C l, W 

Y-axis moment (pitch) 

Cm 

Cm,s 

Cm,wt 

c 

^m,w 

Z-axis moment (yaw) 

Cn 

Cn,s 

^n,wt 



c 

^n,w 


Q 

Subscript P sometimes used to denote principal axes. 




TABLE VI.- SCALE FACTORS a FOR DYNAMIC SCALING 


Area* 5 

Volume 

Velocity ....... 

Acceleration 

Mass 

Mass -flow 

Weight 

Force coefficient . . 
Moment coefficient . 

Wing loading 

Time 

Dynamic pressure . . 

Force 

Angular velocity . . . 
Angular acceleration 

Moment 

Moment of inertia . . 


V 

R 2 
l 

Rg 1 / 2 R ^ 1 / 2 
R g 

r pV 

Rp Rg 1 / 2 R ^ 5 / 2 
Rp Rg R ^ 2 


loO 

1.0 


Rp Rg R^ 

R i 1 / 2 Rg " 1 / 2 

Rp Rg R^ 

Rp Rg Rj 2 
R ^ -1 / 2 Rg 1 / 2 

R^ 1 R g 

Rp Rg R ^ 4 


a 


Scale factor is ratio of model quantity to prototype quantity; for example, 

Model area = R^ 2 x Prototype area 

k Definitions of symbols: 

R _ Model length 
^ Prototype length 



Acceleration of gravity at model altitude 
Acceleration of gravity at prototype altitude 


R 


P 


Air density at model altitude 
Air density at prototype altitude 


107 













Figure 6.- Relation of heading angle H and flight -path angle y to 

earth-centered inertial axes. 
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Figure 7.- Resolution of relative velocity into components along vehicle body axes. 
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Figure 8.- Directions of surface and geostrophic winds. 




z 

Figure 9.- Alinement with respect to vehicle reference axes of thrust and torque 

due to rotating mass of engine. 



Figure 10. - Determination of vertical location of center of gravity. 
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Figure 13.- Determination of vertical location of center of gravity and rolling moments of inertia 

for full-scale airplanes. (Reproduced from ref. 3.) 






Figure 15.- Damped angle -of -attack oscillation assumed in analysis 


Trim angle of 
a ttack 


of appendix D. 


NASA-Langley, 1972 
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